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o ■ 

^ ■ Abstract: We discuss in details a simple, purely bosonic, quantum field theory be- 

^ ■ longing to larger class of models with the following properties: 

■ a) They are asymptotically free, with a dynamically generated mass scale, 
b) They have a space of parameters which gets quantum corrections drastically modify- 
ing the classical singularity structure. The quantum theory can have massless solitons, 

■ Argyres-Douglas-like CFTs, exhibit confinement, etc. . . 

O ■ c) The physics can, to a large extent, be worked out in models with a large number of 

^ . supersymmetries as well as in purely bosonic ones. In the former case, exact BPS mass 



O 



formulas can be derived, brane constructions and embedding in M theory do exist, 
d) The models have an interesting 1/N expansion, and it is possible to define a double 
scaling limit in the sense of the "old" matrix models when approaching the singularities 
in parameter space. 

^ . These properties make these theories very good toy models for four dimensional gauge 

theories with Higgs fields, and provide a framework where the effects of breaking su- 
^ . persymmetry can be explicitly studied. In our model, we work out in details the 

quantum space of parameters. We obtain the non-local lagrangian description of the 
Argyres-Douglas-like CFT, and show that it admits a strongly coupled fixed point. 
We also explicitly demonstrate property d). The possibility of defining such double 
scaling limits was not anticipated on the gauge theory side, and could be of interest to 
understand the gauge theory/string theory correspondence. 

Keywords: Nonperturbative Effects, 1/N Expansion, Sigma Models. 
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1. General presentation 

1.1 Motivations 

The present paper, which extends and provides full details on a previous work [1], 
grew up from the contradictory feelings one might have with regard to the dramatic 
progresses the theory of strongly coupled supersymmetric gauge fields and strings have 
undergone in the past six years (for reviews, see [2], [3], [4]). Though it is likely that 
the general intuition gained in studying supersymmetric examples will be useful to 
eventually understand more realistic theories, the specific methods and the analytic 
solutions of the supersymmetric theories will probably be irrelevant for solving gen- 
uinely non-supersymmetric models. This is due to the fact that those theories cannot 
be viewed as perturbation of supersymmetric theories (at least if the supersymmetric 
theories have at least eight supercharges), as early [5] as well as more recent [6] works 
tend to demonstrate. For example, a gravity approximation to the hypothetical string 
theory description of gauge theories, which can yield useful insights in supersymmetric 
or nearly supersymmetric models, cannot apply to QCD since the hadron spectrum is 
string-like. These limitations have led me to try to find a simple framework where the 
effects of breaking supersymmetry could be analyzed. The present paper is devoted to 
a detailed study of the simplest non-trivial, non-supersymmetric, model belonging to 
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a large class of theories which are simple enough to be tractable even in their strongly 
coupled, non- super symmetric, regime, but complex enough so that many interesting 
questions about gauge theories have a counterpart in the simple models. More precisely, 
and as will be explained in the following, our simple models have all of the following 
general properties: 

a) The models are tractable from the non-supersymmetric versions to the supersym- 
metric ones. In the latter case exact results can be obtained (BPS mass formulas in 
particular) in strict parallel with what is known for four dimensional supersymmetric 
gauge theories. Formulas can actually quantitatively coincide in these cases. Asymp- 
totically free as well as conformal field theories can be studied. 

b) The models have an analogue of a moduli space, with generically both weakly coupled 
and strongly coupled regions. Strong quantum corrections then drastically modify the 
classical structure. At weak coupling we can have solitons playing the role of magnetic 
monopoles or dyons, and these can become masslcss at strong coupling singularities. 
Argyres-Douglas -like CFTs [7] can appear at strong coupling. All these phenomena 
can occur and be studied in both supersymmetric and non-supersymmetric theories. 

c) The supersymmetric versions of the models admit constructions in terms of branes, 
and they can be solved via M theory. 

d) The models have a non-trivial 1/N expansion a la 't Hooft [8]. The large N limit 
can be unconventional, as for = 2 super Yang-Mills [9]. 

We will illustrate a), b) and d) in this article; c) is already known, as we will review 
below. Two additional properties would also be desirable, 

e) The supersymmetric versions of the models can be geometrically engineered as in 
[10]. 

f) The models are dual a la Maldaccna to some kind of string theory. 

Though I am not aware of any explicit construction, it is very likely that e) is true, 
as explained later. As for f), it is plausible that it could be true in view of c) and d), 
but we will have unfortunately nothing to say about this fascinating possibility in this 
paper. 

A fundamental question of principle, that we would also like to address, is whether 
it is possible to prove, or at least to get a good general understanding, that gauge 
theories (supersymmetric or not) or other four dimensional field theories can have a 
description in terms of string theories. The modern starting point is a conjecture 
[11] motivated by the relationship between supersymmetric D-branes and solitons in 
closed string theories. It is not clear whether this intuitive understanding of the gauge 
theory /string theory correspondence makes sense in the general non-supersymmetric 
case. Interestingly, the results of the present paper suggest a way to understand the 
possible proliferation of dualities between four dimensional field theories and string 
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theories. The idea is to show that double scahng hmits [13] can be defined in the 
vicinity of the singularities of the moduli (or parameter) spaces of the gauge theories 
when a non-trivial interacting physics develops at low energies. The double scaled 
theory is then a string theory [12] that can be shown to describe the interacting low 
energy degrees of freedom. We will explicitly demonstrate that such double scaling 
limits can be defined for the model we consider in this paper. These double scaling 
limits have interesting properties, and we will try to provide a more detailed study in 
a forthcoming publication [14]. 

1.2 A family of toy models for gauge theories with Higgs fields 

The fact that good toy models exist for four dimensional gauge theories is of course 
not new. A quarter of a century ago, Polyakov showed in [15] that two dimensional 
non-linear a models can be asymptotically free, undergo dimensional transmutation, 
and develop infrared slavery, the landmarks of interesting four dimensional gauge the- 
ories. Based on this idea, many interesting results were then obtained (see e.g. [16]). 
The models we will propose are only modest extensions of the original non-linear o 
models, the main new input being a way to mimic the presence of Higgs fields, which 
are conspicuous in modern studies. It seems that either it was not known that the stan- 
dard non-linear a models could be modified in a way that would make them suitable 
to compare with modern gauge theory studies, or it was not known that with these 
modifications the models would still be tractable enough in the non-supersymmetric 
cases. It is of course disappointing that we have to restrict ourselves to two-dimensional 
models, but, even at the turn of the millenium, this is still the price to pay to discuss 
the effect of breaking supersymmetry in strongly coupled theories. I hope that the re- 
sults of [1] and of this paper will convince the reader that the models we are proposing 
constitute a nice, if modest, playground to study this fundamental question. 

We now turn to describe the basic idea which led to the construction of our models. 
One important peculiarity of supersymmetric gauge theories a la Seiberg-Witten [17] 
is that they have massless scalar fields and a continuous moduli space of vacua. This is 
heavily used in the discussion of the theories (as well as in many other supersymmetric 
systems), and may appear as being an obvious and impassable obstacle in trying to 
find non-supersymmetric analogues. In non-supersymmetric theories, any vacuum de- 
generacy that may be present classically is generically lifted quantum mechanically. In 
two dimensions, the situation is even worse: even supersymmetric theories cannot have 
a continuous moduli space, because the strong infrared fiuctuations always make the 
vacuum wave functional to spread over the whole would-be moduli space (a discrete 
set may remain [18]). However, it turns out that this is only an outward problem, for 
the following reason. In four dimensions, the moduli space is parametrized by Higgs 
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vacuum expectation values (or more precisely by the vevs of independent gauge invari- 
ant combinations of Higgs fields) . The physics is interesting because the masses of the 
gauge bosons, which govern the low energy coupling, depend on the Higgs vevs via 
the usual Higgs mechanism. Moving on the moduli space is then equivalent to varying 
the low energy coupling, and a very interesting physics is associated to the transition 
from weak coupling to strong coupling: appearance of strong coupling singularities [17], 
non-trivial CFTs [7], rearrangement of the spectrum of stable states [19], etc. . . All this 
physics is largely independent of the fact that the Higgs potential has fiat directions, 
but strongly depends on our ability to vary the low energy coupling. We are thus led to 
the conclusion that a good non-supersymmetric analogue of the moduli space of super- 
symmetric theories could be a space of parameters on which the low energy coupling 
depends. This is too vague, since there are a priori many ways to vary the low energy 
coupling, for example by putting the theory on a sphere of varying radius or considering 
a finite temperature. To stick as close as possible to the supersymmetric gauge theory 
case, we will in general consider parameters that correspond to giving masses to the 
fields that contributes with a minus sign to the f3 function. In the case of the non-linear 
a model with target space the — 1 sphere S^~^ that we will consider in the present 
paper, this simply amounts to giving a mass to the N — 1 would-be Goldstone bosons. 
These mass parameters play the role of Higgs vevs, and span a space of parameters 
A4. One of our main goal is to compute the quantum corrections to M., in the same 
sense as quantum corrections to the moduli space of supersymmetric gauge theories 
were computed in [17] or [20]. 

Though the very simple idea presented above certainly suggests that the space 
of mass parameters in non-linear a models is an interesting object to consider, the 
reader may not be convinced that it is really a good analogue of the moduli spaces of 
supersymmetric gauge theories. The present author himself actually became convinced 
that it is the case only after he became aware of papers by Hanany and Hori [21] and 
Dorey and collaborators [22], where the M — 2 supersymmetric £,P^~^ non-hnear a 
model with mass terms is discussed (in this context, the mass terms are called "twisted 
masses," and are in one to one correspondence with the holomorphic isometrics of the 
target Kahler manifold [23]). In [21], a brane construction of this supersymmetric model 
is given, and M theory is used to obtain non-perturbative results, in strict parallel to 
what was done in the case of A/" = 2 super Yang-Mills [24]. The striking, in some sense 
quantitative, similarity with M = 2 super Yang- Mills was then further discussed in [22]. 
The property making brane constructions possible is that the mass parameters can be 
interpreted in the <CP^~^ model as vector multiplets vevs in a gauged linear formulation 
of the model, and thus literally parametrize a Coulomb branch. The works [21, 22] thus 
clearly demonstrate that an excellent toy model ior J\f — 2 super Yang-Mills in four 
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dimensions is the N — 2 <CP^~^ non-linear a model with mass terms in two dimensions 
(in spite of the fact that jV" = 2 in four dimensions corresponds to eight supercharges, 
whereas jV = 2 in two dimensions corresponds to four supercharges). In retrospect, 
this analogy is not too surprising, and there are actually many independent arguments 
in favor of this correspondence. For example, both types of theories are known to have 
similar non-renormalization theorems ([25], [26]), and both admit topological twists 
([27], [28]). 1 think that this latter property is particularly significant, because the exact 
results a la Seiberg-Witten ([17], [2]) arc likely to have a scmi-topological origin, very 
much like the exact results one can obtain in two dimensions [29]. Another important 
fact is that effective superpotentials in two dimensions are very similar to effective 
prepotentials in four dimensions. This was emphasized for example in [30] where two 
dimensional superpotentials were geometrically engineered using singular Calabi-Yau 
fourfolds and mirror symmetry, in the same spirit as four dimensional prepotentials can 
be engineered using singular Calabi-Yau threefolds and mirror symmetry [10]. Though 
the models considered in [30] were different from the one we are proposing in this paper, 
a similar construction is probably possible. 

We believe that this analogy between supersymmetric theories in four and two 
dimensions can be fruitfully extended to the non-supersymmetric cases, and used to 
study the effects of breaking supersymmetry. We present below a short dictionary for 
the correspondence between four dimensional and two dimensional models. Some of 
the entries will be exemplified in later sections. 



4 D gauge theories 


2 D non-linear a models 


Gauge bosons 


Goldstonc bosons 


Gauge group 


Isometry group 


Number of colors 


Dimension of the compact target space 


Gauge coupling constant 


Inverse radius of the target manifold 


Higgs fields 


Mass terms 


Moduli space M. 


Space of mass parameters Ai 


Monopoles and dyons 


Kinks 


S duality 


Kramers- Wannier duality 


Argyrcs-Douglas CFT 


Ginzburg-Landau CFT 


Eight supercharges 


Four supercharges 


Effective prepotential T 


Effective superpotential W 


String description 


Branched polymer description 



Table 1: A gauge theory/cr model dictionary 
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1.3 Plan of the paper 



We give a general discussion of mass terms, including their rcnormalization properties, 
in Section 2.1. We work out a few simple examples when the target space is the sphere 
S^^^, including a model that we will call after C. Neumann and on which we will 
focus. We also briefly discuss the supersymmctric generalizations. In Section 2.2, the 
semiclassical properties of the Neumann model are explored, including the singularity 
structure of the classical space of parameters M.ch bound states and solitons. Section 
3 is devoted to the expansion. Emphasis is put on the peculiarities introduced by 
the fact that the dimension of A4 is of order N, and on general limitations of the 1/N 
expansion in our model. Wc compute the quantum corrections to the metric as well 
as the mass of some stable particles. We then show how to go beyond the standard 
expansion, and wc obtain a provisional picture of the singularity structure of the 
quantum space of parameters Aiq (also simply denoted by Ai in the following). In 
Section 4 we work out the full structure of in the large N limit. In Section 5 we 
analyze the theory in the vicinity of singularities on M., for any finite N. We show 
that the non-trivial infrared physics is most naturally described in terms of a non-local 
lagrangian analogous to a theory of electric and magnetic charges. This brings the 
analogy with an Argy res- Douglas theory to a climax. Due to the simplicity of two 
dimensional theories, wc are able to show that the non-local theory has a strongly 
coupled fixed point, and we recover the results of Sections 3 and 4 independently of a 
large N approximation. In Section 6, we briefly show that double scaling limits, in the 
sense of [13], can be deflned as the singularities on Ai arc approached. This has never 
been studied in the context of two dimensional non-linear a models, and our results 
suggest very interesting possibilities for four dimensional gauge theories that have not 
been anticipated. Moreover, these double scaling limits have interesting properties, 
which we will discuss elsewhere [14] . 

We have also included three appendices. In Appendix A, we find the most general 
static finite energy solutions to the field equations of our model at arbitrary A^. A 
notable result is that wc obtain static solutions describing several standard sine-Gordon 
solitons at arbitrary distances from each other. In Appendix B we derive several simple 
formulas used in the main text to study the large limit of the theory. Finally in 
Appendix C we compute the 1/A^ corrections to the equation of the critical hypersurface 
obtained in Section 3. The calculation illustrates some generic properties of the 1/A^ 
expansion, as well as some subtleties associated with infrared divergences that can 
plague this expansion in our model. The result gives a consistency check of the simple 
calculations of Section 3. 

To the opposite of [1], we provide in this work detailed derivations and elementary 
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discussions of various points. This explains in part the length of the paper. 



2. The Neumann model 

2.1 Mass terms in the non-linecir sigma model 

In the following, the target space of our non-hnear a models is the N — 1 sphere S^~^. 
The discussion could be easily generalized to sigma models on symmetric spaces, for 
example. At the technical level, the renormalization properties of the two dimensional 
non-linear sigma model are non-trivial [31], and their renormalizability can be proven 
with the help of a sort of Zinn- Justin equation, very much like the renormalizability of 
gauge theories is analyzed (see e.g. [32]). In the latter case, it is not straightforward 
to add mass terms for the gauge bosons that do not spoil the renormalizability of the 
model. To do so, it is necessary to introduce scalar fields, the Higgs bosons, and the 
gauge bosons masses are then largely determined by the transformation properties of 
the Higgs fields under the gauge group G through the Higgs mechanism. The situ- 
ation is similar for mass terms in non-linear sigma models: they generally spoil the 
renormalizability of the theory, and from this point of view they are most naturally 
characterized by their transformation properties under the isometry group G = 0{N) 
of the target space. The explicit form of the mass terms in the lagrangian is then 
determined unambiguously by G invariance and power counting. 

2.1.1 Renormalization theory 

Working in the euclidean, the lagrangian without mass terms is 

I iv-i 

i-kin = 2 9ij da^ida^j, (2.1) 

where gij is the standard 0{N) invariant metric on a sphere of radius g being the 
dimensionless coupling constant analogous to the gauge coupling constant. By taking 
the $iS to be the usual cartesian coordinates, we have 

9i3 = + , ^TSU^, (2-2) 

The model is defined by a path integral over the fields $i, l<i<A^ — 1, with an 
O(A^) invariant measure. We have interactions of the form 51^" (9$)^ for all n > 1. 
A mass term is defined to be a couple (F, d) where T is an irreducible representation of 
SO(A^) and d is the canonical dimension of the mass parameters. In two dimensions, 
it is consistent to have d — 1 or d — 2. In the latter case, the parameters actually 
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correspond to mass squared. Different choices of T and d do not necessarily correspond 
to independent mass terms, as we will see. li d — 1 the mass parameters are denoted 
by m = mi^...i^, and if d = 2 by h = hi^...i^. The ijs are SO(A^) indices, and the tensors 
m or h satisfies some constraints depending on T. 

Instead of working with the lagrangian (2.1), which is possible but awkward, we 
will introduce a Lagrange multiplier field a and work in a representation where the 
0{N) symmetry is linearly realized: 



^kin 



+ ('$2_ 1 , ^2.3) 
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with $ = ($1, . . . , $jv). Ehminating from (2.3) by using the constraint enforced 
by a, one recovers the non-linear lagrangian (2.1). We will regulate the theory by a 
simple momentum cutoff Aq, which is manifestly 0(iV) invariant. Introducing m or h 
into the game, the action must then be taken to be the most general relativistic local 
functional invariant under O(A^) (by varying both the fields and the parameters) and 
compatible with power counting. The canonical dimensions of the fields are 

[$] = 0, H=2. (2.4) 

As we will further discussed in Section 3.1, these dimensions are relevant to study the 
UV finiteness of both ordinary perturbation theory and of the 1/A^ expansion because 
the model is asymptotically free. Due to its canonical dimension, a can only appear in 
a term of the form F(#^) ct, where F is an arbitrary function. By integrating over a, 
we get the quantum constraint 

F($2) = 0, (2.5) 
which can always be solved in the UV as 

*^ = 4?' 

where ZZg is a renormalization constant. This shows that, as long as we integrate over 
a without introducing a source for this field, F can be taken without loss of generality to 
be proportional to — l/{ZZgg'^), and can be replaced by its constant value in the 
other terms of the lagrangian. The subtleties associated with the renormalization of a- 
dependent quantities will only show up in Appendix C, and is discussed there. Without 
any mass term, the only non-zero dimension two operator that remains is proportional 
to da^da^, and from this we deduce the form of the renormalized lagrangian [31] 

Lkin = - da^da^ + - « ( ^ | ■ (2.7) 

2 2 V ZZgQ^ J ^ ' 
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$ and g are now the renormalized fields and coupling constant, related to the bare 
quantities by 

In the following, we investigate the simplest mass terms and find the scalar potential 
V they correspond to. 

Singlet Because of the constraint (2.6), singlet mass terms are trivial and simply 
correspond to adding a constant to the lagrangian. 

Vector hi of dimension two This term corresponds to a standard magnetic field, 

^(.,2) = -^(.,2) h$ = -Z(,,2) /i$iV, (2.9) 

where we have used the O(A^) symmetry to ahgn the magnetic field with the A'"th 
direction. This is a linear symmetry breaking term, and thus ^(^,2) — 1- By eliminating 
we have 



V{v,2) = -Z(v,2) ^ 



7V-1 



1-^2 ^$2. (2.10) 



i=l 



In addition to giving a mass 2)fi'^ to the would-be Goldstone bosons, this term 

also produces an infinite sum of new interactions. This model can easily be studied in 
the large N limit, but it is not particularly interesting. 

Vector rrii of dimension one The only new invariant operators of dimension 2 are 
(m^)^ and m^. One has also the dimension one operator m^. The potential is then 

V^v,i) = (m*)2 - Mm* - Z(;,)m2 

= „,2^^ _ z'^^^^^ ^^^^ _ ^^^^2_ (2.11) 



M is a new mass parameter that generically enters the problem. However, it can be 
consistently set to zero thanks to the symmetry $jv ^ —^n- ^[vi) renormalizes the 
vacuum energy. 

Antisymmetric tensor of dimension two No new operator can be constructed, 

V^a,2) = 0. (2.12) 



10 



Traceless symmetric tensor hij of dimension two This corresponds to 

VW, = -%^*h# = -%2lf:W?. (2.13) 

i=l 

By diagonalizing the matrix hij, we see that we have N — 1 independent parameters, 
since the hiS satisfy X^ili — 0- This tracelessness condition can actually be waived 
without changing the physics, since the trace part of h corresponds to adding a constant 
to the lagrangian. This model, that we will call the Neumann model for reasons to 
become clear later in this Section, is the one on which we want to focus in this paper. 
We will usually use the — 1 independent variables 

Vi^hN-hi, l<i<N -1. (2.14) 

Note that in the particular case vi — • • • — vn-i — rn? we recover the model for a 
dimension one vector mass parameter with M = 0. By eliminating $jv we find, up to 
a constant term, the very simple potential 

ry N-l 

V = ^Y.->^1 (2.15) 

i=l 

Antisymmetric tensor aij of dimension one Dimension two operators can be 
constructed by using tijki = dijdki- This tensor decomposes into four irreducible repre- 
sentations of SO(A^), amongst which only two can be used to construct invariants with 
the help of the $,jS: the trivial representation tra^, and tij-j}-. We thus get a special 
case of the previous model, for which 

N 

hij = - ^ aikajk- (2-16) 

This special case is particularly significant, however, because it is on this form that the 
Neumann model can be supersymmetrized [33]. 

The Neumann model is singled out by the fact that it is the most general model 
with only quadratic terms in the fields $ when the Lagrange multiplier a is introduced. 
This makes the analysis of the large limit very simple as we will see in Section 3. More 
general models can nevertheless be studied in the large limit, with the ideas presented 
in the present paper, but at the expense of introducing additional auxiliary fields (the 
large N limit can be subtle in some cases, though). In this sense the Neumann model 
is "minimal," since it can be studied by introducing the minimal number of auxiliary 
fields. 

It is natural to consider another simple model. 
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Traceless symmetric tensor rriij of dimension 1 To work out the potential term, 
we must decompose tijki — mijmki into irreducible representations of SO(A'"). Four 
such representations occur, three of which can be used to construct three independent 
invariants, including the trivial representation which renormalizes the vacuum energy. 
This implies that, in addition to the mass parameters themselves, the model depends on 
a new dimensionless coupling constant G. Generically, we will also have the dimension 
one operator $m$ and a new mass parameter M. The potential is 

V^^^-y^ = ^fiil ^m^* — (#m#)2 - 1) M$m$ - Z'(^^y^ trm^ 

i=l \i=l 

TV TV 

-^(M) E - ^(M) E <■ (2-17) 

M can be set to zero thanks to the symmetry i— > — m^. G can also be consistently 
set to zero, since we then recover a special case of the Neumann model, for which 

TV 

hij = - ^ miknijk. (2.18) 

k=l 

Considering G 7^ is however interesting, since the model can be supersymmetrized 

for G = g [33] (in the non-supersymmetric version, it is not consistent to set G = g, 
and G and g are then two independent coupling constants, see below). We thus have 
two natural supersymmetric versions of the Neumann model, with h taking the special 
form (2.16), or (2.18) with G = g. In the latter case, the scalar potential can be written 
most elegantly by introducing a new auxiliary field a, which turns out to be in the same 
supersymmetry multiplet as the Lagrange multiplier a, 

V,.,i, = ^ + (2-19) 

1=1 

This is the M — 1 supersymmetric version of the M — 2 potential considered in [21] 
and [22], which is 

1 ^ 

where a and the $iS are now complex fields. The Af — 2 theory with the potential 
(2.20) shows quantitative similarities with Af — 2 super Yang- Mills in four dimensions 
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[21, 22], as we have already pointed out in Section 1. We thus see that the Neumann 
model can be viewed as a bosonic version of the supersymmetric theories studied in 
[21, 22]. 

2.1.2 Renormalization constants 

One-loop formulas for the various renormalization constants can be easily obtained 
with the background field method and using Ricmann normal coordinates, which is 
very elementary in the case of the sphere. We recover the well known formulas 

Z, = l + ^/ln-^ + 0(^^), (2.21) 
^ = l + |^ln£ + 0(^'), (2.22) 



and we also obtain 

%2) = l-|^ln^ + 0(^?^). (2.23) 

Aq is the UV cutoff and fi an arbitrary renormalization scale. The RG functions of 
the Neumann model are then (we have indicated /3 up to two loops [34], since this is 
needed later) 

The theory is asymptotically free for > 3 [15], with a dynamically generated mass 
scale A defined at one loop by the equation for the couphng at scale 

1 N -2 a ^^^^^ 
In^- (2.27) 



g^ili) 2n A 

We have also computed 



^2 _|_ 4(^2 

Zis,i) = 1 - In ^ + Oig% (2.28) 

= + (2.29) 



13 



from which we can deduce the /? function for the couphng G, 

PGiG^g')^- ^'^^^'-^''K oig% (2.30) 

71 

This shows unambiguously that C is a new couphng, independent of g, and thus that 
introducing mass terms for the Goldstone bosons will generically introduce new dimen- 
sionless coupling constants into the theory. This is another common point with gauge 
theories, where in addition to the gauge coupling constant we have the couplings in the 

Higgs potential. However, in two dimensions, the new couplings cannot spoil asymp- 
totic freedom: in spite of possible plus signs in f3 functions like (2.30), the physical 
coupling is really something like G^w? where m is some mass parameter, and is always 
irrelevant in the UV. 

2.2 The Neumann model 

Prom now on, we focus exclusively on the Neumann model, whose minkowskian classical 
lagrangian and field equations are 

{^^^^'^^ + - f - ^ j ' (2-31) 

d^d^^, = {h, - a) (2.32) 
The space of parameters A4 is spanned by the N — 1 real independent variables 

Vi — hjsr — hi. (2.33) 

A singularity on is a point where some of the degrees of freedom become massless. 
The set of all the singularities on M. is called the critical hypersurface H C M.. We 
will distinguish between the classical critical hypersurface Hci and the quantum critical 
hypersurface Tiq. We will also speak loosely of a classical M.ci and quantum M.q space 
of parameters, meaning M. equipped with Hd or Hg respectively. For generic values of 
the parameters Vi, the model has a symmetry 1^2(1) x • • • x Z2(iv) with 

Z2w:^,-^(-l)'^^-^,-- (2.34) 

When the ViS coincide for p disjoint subsets of indices Ji, . . . , /p, of respective cardinality 
ki, . . . , kp, the symmetry group is enlarged to Xj^/j(j...|j 7^22(1) x 0{ki) x • • • x 0{kp). 

The Neumann model, in addition of being a natural extension of the standard non- 
linear a model, is also a generalization of the sine-Gordon model. Indeed the classical 
potential 

1 ^ 

Ki = -^E^^^' (2.35) 
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reduces for = 2 to the sine-Gordon potential 



VsG^ 



4^2 



COS 26 



(2.36) 



by writing 



$1 = 



sin^ 



$2 = 



COS^ 



(2.37) 



9 



9 



In the case = 3, our model is also related, through Haldane's map [35], to the 
anisotropic XYZ quantum large integer spin chain in one dimension, or equivalently 
to the anisotropic classical Heisenberg model in two dimensions. This analogy is quite 
helpful to understand the physics of the model. The model for = 3 has instantons, 
and we could introduce in this case a non-zero 9 angle. The physics would then strongly 
depend on 6. We will not do that in the following, however, since we want to study 
the models for arbitrary in a unified way. Considering a non-zero 6 would be more 
natural in the context of a work on the CP^ model with mass terms. 

The model defined by (2.31) is also the field theoretic generahzation of a famous 
integrable system in classical mechanics corresponding to the motion of a particle on a 
sphere with a quadratic potential. 



Prom This problem was first studied by C. Neumann 150 years ago (for N — 3), hence 
the name "Neumann" for the model (2.31). The time independent solutions to the 
field theory satisfy (2.38). In Appendix A we use integrability to find the most general 
solitonic solutions of the Neumann model, generalizing the famous sine-Gordon solitons. 

2.3 The weakly coupled theory 

At weak coupling, the classical vacua of the theory correspond to the minima of the 
potential (2.35) defined on the sphere Xlili ^1 ~ Vfi'^- restrict ourselves, without 

loss of generality, to the region A4n of parameter space defined by Vi > 0, then we have 
generically two equivalent vacua related by the spontaneously broken 7i2(N) symmetry, 
such that 



By expanding around one of these vacua, we read from the lagrangian that the A^ — 1 
"Goldstone bosons" $i, l<i<A^ — 1, have masses 




1=1 



(2.38) 




(2.39) 



rrii = ^/v'i. 



(2.40) 
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This relation together with (2.27) shows that the coupUng is weak in the region 



t;^ > for 1 < i < - 1. (2.41) 

It is actually sufficient to have A^— 2 heavy fields for the physics to be weakly interacting, 
for example 

> A^ for l<i<N -2. (2.42) 

The low energy theory is then a sine-Gordon model whose elementary field has a mass 
y/VN-i and whose non-running coupling is small. 

The classical hypersurface of singularities Hci in Al at is trivially determined by 
(2.40) to be the union of the hyperplanes Vi — 0. In addition to A4n, there are A" — 1 
other regions A4j C A4 defined by 

Mj = {{vi,..., VN-i) I vi'^ = Vi- Vj > 0}, (2.43) 

or equivalently by hj = max, hi in A4j. Hd H A4j is then the union of the hyperplanes 
v^^^ = in M.j. We have reprensented Aid for the cases N — 3 and A" = 4 in Figure 1. 




Figure 1: The classical space of parameters Md in the cases N = 3 (left) and N = 4 (right). 
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2.3.1 Integrability? 

In view of the relation of our model with the Neumann integrable system, the O(A^) 
non-linear a model and the sine-Gordon model, the reader may wonder whether it 
could be integrable itself or not. However, the simple arguments that can be used 
to prove integrability when v\ — ■ ■ ■ — I'^v-i = or A?" = 2 no longer work, because 
the model combines in general the complications due to both an explicit mass term 
and a non-zero /9 function. In other words, there are enough operators to spoil the 
higher conservation laws of the standard non-linear a model. Integrability can actually 
be easily disproved in the weakly coupled region of parameter space, for example by 
computing the S matrix element for the scattering of elementary particles, 

{h qi) + (i, q2) {k, pi) + {I, P2), (2.44) 

where qi = (cUj^gj, gi), q2 and Pi,P2 are the incoming and outgoing two-momenta re- 
spectively. In the Born approximation, which is valid at weak coupling, we have 

{k,pi;l,p2 \S-1\ i,qi;j,q2) = (^^^Hpi + P2 - Qi - q2) 

[s 6ij6ki + 1 6ikSji + u 6ii6jkj (2.45) 



where s = (qi + q2)^, t — {qj — pi)^ and u — {qi — p2)^ are the usual Mandelstam 
variables, and with the normalization 

{j,P \hQ) = uJi,qSijS{q-p) = \lq^ + mj5ij5{q-p). (2.46) 

It is clear from (2.45) that processes where the number of particles of a given mass 
changes are allowed, which proves that the model cannot be integrable in this regime 
[37]. In the 0(A^ — 1) symmetric case Vi = ■ ■ ■ = v^^i where all the elementary particles 
have the same mass, the S matrix element can be written 

{k,pi;l,p2 1-5-11 i,qi;j,q2) = u;q,5(qi - pi) u;q^5{q2 - P2) 

{SijSki Si + SikSji S2 + SiiSjk S3) + {k l,pi P2), (2.47) 



with 

On this form, it is clear that the Yang-Baxter equation is violated [37], and thus the 
theory cannot be integrable even in the most symmetric case. The above reasoning 
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does not exclude that the model could be integrable for some special values of the 
parameters at strong coupling, but apart from the case vi — ■ ■ • — vn-i = we think 
it is very unlikely. This is confirmed by the analysis in the large N approximation, see 
Section 3. 

2.3.2 Bound states 

From what is known for the sinc-Gordon model, we can suspect that we have bound 
states of the elementary particles at weak coupling. This is a priori possible because the 
interactions between the $jS are attractive thanks to the derivatives in the interaction 
terms (see (2.1) and (2.2)). At weak coupling, we can use a non-relativistic approxi- 
mation to investigate the bound state spectrum. The quantum mechanical attractive 
potential in the two-particle subspace can be straightforwardly deduced from the S 
matrix element (2.45) and is simply 

V{X)^~g^S{X)^J (2.49) 

where J acts in internal space, 

{M\J\ ij) = SijSki. (2.50) 

We see that processes like (ii) — > (jj) with i j are possible, as long as the kinematical 
non-relativistic constraint = nij is satisfied. In general, stable bound states will exist 
between particles of the lowest mass. Let us study the two-particle bound states in the 
symmetric case Vi = ■ ■ ■ = v^-i = v. This amounts to solving the Shrodinger equation 
for two particles of mass a/u interacting through (2.49). The wavcfunctions ipijlxi, X2) 
must satisfy ipalxi, X2) = 'ipii{x2, Xi) due to bose statistics. It is very elementary to see 
that there is a unique bound state which is a singlet of 0{N — 1) {ipij oc 5ij) and whose 
mass is 

mb = 2V^(l--(iV-l)V). (2.51) 

The symmetric mixing between the — 1 two-particle states making up the 

bound state dramatically stabilizes it. This is particularly significant for the physics of 
the model, as we will see later. In particular, in the large N limit a la 't Hooft [8] 

N ^00, g ^0, g^N = constant, (2.52) 

the interactions between the fields $i are of order 1/N, but they can nevertheless form 
a bound state whose binding energy is of order N^, as (2.51) shows. We will discuss 
the properties of this bound state further in Section 3. It is also natural to look for 
multiparticle bound states, which we know must exist in the limit where the model 
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approach the sine-Gordon model (when one of the ViS is much smaller than the others). 
This is non-trivial, because the J matrices corresponding to different pairs of particles 
don't always commute, and we have not tried to study the corresponding Schrodinger 
equation. Note that these multiparticle bound states will be unstable in the field theory, 
which is unlike the case of the sine-Gordon model, but their lifetime can be large. In 
the case where the ViS are not all equal, we will still have a stable bound state in which 
the pair $-$ corresponding to the particle of lowest mass dominate. We will see how 
(2.51) generalizes in this case in Section 3 in the large N limit. 

2.3.3 Solitons 

The most general time independent, finite energy solutions to the field equations (2.32) 
are derived in Appendix A. They satisfy (2.38) and must tend towards one of the vacua 
(2.39) of our theory. Let us assume for the moment that Vi > V2 > ■ • •i'at-i- The 
only solutions corresponding to stable particles are then found by restricting the fields 
to $1 = ■ ■ ■ = $Ar_2 — 0, and correspond to the standard sine-Gordon soliton and 
anti-soliton of mass 

rrisoi = — (2.53) 

9 

For example the soliton solution joining the North pole of the sphere S^~^ at x = — oo 
to the South pole at x — +oo and centered at x — Xc is given by 

sin 9 ^ cos 9 ^ ,s 

^N-i = ' = ' (2-54) 

9 9 

with 

tan^ = eV^(^-^-). (2.55) 
There are also solutions corresponding to sine-Gordon solitons of masses 



rriiM = — (2.56) 
9 

when the only varying fields are taken to be $jv and The particles associated with 
these solutions are unstable and would decay to the stable solution (2.54), (2.55) by 
emitting elementary quanta. 

There are also more general solutions, which are worked out in details in Appendix 
A, describing a succession of sine-Gordon kinks and anti- kinks of different masses. The 
force between these kinks magically cancel at the classical level, a consequence of the 
integrability of the equations (2.38). It is very unlikely that this property is maintained 
at the quantum level. The static solution for the case N — 3 is depicted in Figure 2. 
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Figure 2: A static solution for the case N = 3, vi = 1.5, V2 = 1, and g = 1, describing 
two kinks of respective masses 2yjv\jg^ and 2y^/(7^, at different values of the distance 
Axc = Xc2 — Xc\ between their centers. The energy density p is represented in thick solid 
line, and the spherical angles and (/> in thin solid line and dashed line respectively. The 
maximum energy density of a given isolated kink is Vi/g^. The formulas used to obtain the 
figure are given in Appendix A. 3. 



When several of the ViS coincide, the soliton solutions have collective coordinates 
correponding to the enhanced 0(p) symmetries. Let us discuss briefly the most sym- 
metric case f 1 = ■ ■ ■ = vn-i = v. The general solution $j^soi(a; — Xc; ^a) is obtained from 
(2.54) and (2.55) by applying an arbitrary SO(A^ — 1) rotation R{^a)- The angles can 
be taken to parametrize the coset SO(A^ — 1)/S0(A^ — 2) = S^^^. The quantization in 
the one-soliton sector, in the moduli space (low energy) approximation, then proceeds 
in the usual way, by restricting the fields to be of the form 

$,(x,t) = <l>,,soi(x-x,(t);rW), (2.57) 

and replacing this ansatz in the lagrangian (2.31) to find the quantum mechanics gov- 
erning the collective coordinates S,'^{t) and Xc{t). The resulting lagrangian turns out to 
be ^ 

^moduli = -"^sol + 2 "^sol xl + Qo^p C^'^, (2.58) 
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where gap is the 0(N — 1) invariant metric on the sphere S^~^ of radius unity. The 
Schrodinger equation corresponding to (2.58) involves the Beltrami laplacian on the 
sphere S^~^, whose eigenvalues and eigenvectors are well-known. We get a tower of 
states of mass 

m^o\,j = rn^oi + J{J + N - 3) (2.59) 

fiUing multiplets of SO(iV — 1) corresponding to the completely symmetric and traceless 
tensors of rank J. The low-lying states must be stable at weak coupling, since then 
v/rusoi — \fv g'^ /2 <C \/v, the mass of the elementary particles. 

We also have soliton/anti-soliton bound states described by the breather solution 
of the sine-Gordon model. It could be interesting to investigate, for example in the 
case f 1 = • • • = fAT-i, the possible relationship between these two-solitons bound states 
and the multiparticles bound states formed by the $iS. We have not performed this 
analysis, however; for our purposes, only the stable $-$ singlet will be relevant. 

3. The 1/N expansion 

3.1 General properties and limitations of the Icirge N approximation 

The large N expansion a la 't Hooft (2.52) is one of the main non-perturbative tool at 
our disposal. The basic idea (for reviews see for example [38]) is to integrate exactly 
over the fields $i in (2.31), which yields an effective action proportional to N . The 
1/A^ expansion is then nothing but a loop expansion for this effective action, and it is 
non-perturbative with respect to the other parameters of the theory. For example, the 
connected vacuum amplitude can be written in the Neumann model as 

00 

W = Y,N^~''^i{'^il^)- (3-1) 

1=0 

A 2n-point function would have an additional factor of 1/N'^~^. Though it is very 
useful, the expansion has nevertheless two important limitations that will bother 
us in the following, and that are likely to show up in gauge theories as well. 

The first limitation comes from the fact that our space of parameters 7W is — 1 
dimensional. In the large limit, instead of working with a discrete set of dimensionless 
parameters Xi = Vi/A, one should really introduce a density 

N-l 

1=1 

and consider only smooth enough functions r{x). What this means in practice is that 
we cannot study in this approximation scheme non-generic parameters where a small 
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number of masses are adjusted to take particular values. To be concrete, let us suppose 
that we want to investigate the region in parameter space where one of the ViS, say vi, 
goes to zero. In the expansion, the ViS typically appear through combinations of 
the form 



We see that, though V is generically of order iV°, it is of order N when vi — > 0. This 
means that the 1/N counting is completely modified in this limit, and the standard 
expansion can no longer be trusted. This problem is of course a serious drawback, 
since we know in particular from gauge theory studies that interesting phenomena 
are usually associated with fine tuning of parameters ("critical points"). Fortunately, 
there are two ways out of this problem. First, one can consider only p dimensional 
subspaces Sp C A4, with p <^ N, for example by choosing Vi = ■ ■ ■ W[{Ar-i)/p] = Wi, 
f [(Ar_i)/p]+i = • • • = V[2{N-i)/p] = W2, etc. . . Thc physlcs on Sp, including special points, 
can then be studied reliably in thc large N expansion. This trick is useful in the 
supersymmetric generalizations of the Neumann model [33]. In the present case, we 
need to do better, however, because interesting physics can be associated with only 
one parameter going to a special value, like Vi — > 0. The idea is then to treat exactly 
the small number of degrees of freedom that play a special role in the particular limit 
we are interested in. For example, in the limit vi — > with all the other f jS positive, 
we will integrate over the N — 2 fields $2, ■ ■ ■ , ^at-i in (2.31), while keeping explicitly 
the fields $1 and <I>jv together with a. More generally, it is actually very convenient 
to always keep explicitly the field $jv when working on Ai^ = {vi > 0} for example, 
since this field is singled out by the fact that {^n) 7^ at weak coupling in this region. 

The second limitation of the 1/N expansion that we will encounter comes from the 
fact that the terms in the expansion can behave badly for some values of the parameters. 
For example, it could be that some of the coefficients Wi{vi/A) in (3.1) are singular 
for some (generic) values of the ViS. Typically, these coefficients can be expanded 
in powers of the coupling constant at weak coupling, but when the coupling grows 
the series might become singular. Beyond that point, the coefficients might again be 
expandable in terms of some "dual" coupling. At the singular points, the standard 1/A^ 
expansion fails completely. We will see in Section 3.3 that the physics can nevertheless 
be extracted in the large N limit, by resumming in some sense the 1/N series. We will 
also see in Section 6 that in our model the singular behaviour of the coefficients of the 
1/N expansion is very specific, and that this has some very interesting consequences. 
The same kind of phenomena are likely to occur in gauge theories as well, for example 
in the vicinity of an Argyres- Douglas point. 




N-l 



(3.3) 
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Renormalization At the technical level, the 1/N expansion needs to be renormal- 
ized. It is actually easy to find the divergent part of the renormalization constants 
that make the theory finite to all orders in 1/N. This is possible because the theory 
is asymptotically free, and the expansion is non-perturbative in g. The divergent 
contributions can then be obtained by summing up ordinary perturbation theory by us- 
ing the renormaUzation group. As is well known, these divergent terms are completely 
determined by the two-loop (3 function and the one-loop 7 and a functions, all given 
in (2.24), (2.25) and (2.26). More precisely, we have the exact formulas 

— — m — H mm — -I- finite terms, (3.4 

g'g 27r A 27r A ^ ' 

\nZ — —— — - Inln — ^ -|- finite terms, (3.5) 
lnZ^s,2) = j;^ — - Inln-^ -|- finite terms, (3.6) 



which yield 

+ 0{1/N^), (3.7) 
Z(s,2) = l + ^lnln^ + ^ + 0{1/N'), (3.8) 

where k and k' arc finite constants. The finite terms in In Z and In ^(^,2) will contribute 
to the divergent terms in Z and Z(^g^2) at order and higher, however. 

3.2 Mcissless states 

To keep the formulas as simple as possible, let us begin by restricting ourselves to the 
0{N — 1) symmetric case vi — ■ ■ ■ — vn-i — v >0. We can then take without loss of 
generality — v and /ij = for 1 < i < — 1. By integrating over the N — 1 fields 
1 < ^ < A^ — 1, in (2.31), and Wick rotating to the euchdean, we obtain, in the 
large A^ limit, the non-local effective action 

S,s[a,^N]^ Ns,s[a,ip] (3.9) 

with 



and 
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^N = VNip (3.10) 

f 72 f ^ a — V 2 a \ 
Seff[a, (p\= J d xi- da^pda^p H — ip - 2^ 1 + sH- (3.11) 
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The functional s[a\ is defined and studied in details in Appendix B, to which we refer the 
reader, and g is the renormahzed couphng (2.27). The effective potential corresponding 
to Seff is 



a — V 2 a a /o -.r.\ 

Ves{a, cp) = ^- cp - ^ In (3.12) 

and the saddle point equations governing the physics in the limit N —>■ oo are 

{a,-v)^, = 

2 1 (3.13) 

At weak coupling v ^ A^, the 'Z'2(n) symmetry is spontaneously broken and ip^ ^ 
(2.39). We will loosely call the "weak coupling" region in parameter space the whole 
region where Z2(jv) is broken and where we have 

{a* = V 
9 1 , (3-14) 

^* ^ i^^^^A^- 

Equations (3.14) give 

{^%) = ^ ^, = 4-' (3-15) 
g^ifx = V^) glff 

which is simply the classical formula (2.39) where the bare coupling has been replaced 
by the effective low energy coupling ges, the running coupling (2.27) evaluated at the 
scale ^yv corresponding to the mass of the elementary particles. On the other hand, 
when V = 0, we know that ($jv) = as a consequence of the SO(A^) symmetry. Thus, 
there must be a phase transition at strong coupling for some f = fc ~ A^, where 
the Z2(Ar) is restored. We would like to understand the nature of this transition, and 
whether i>c = or > 0. Since it occurs at strong coupling in a non-abelian system, 
the standard intuition would suggest that a mass gap is formed and thus that any phase 
transition must be first order. This would be the case for example if Vc — 0. However, 
we are going to show that this is not what happens: the transition is second order, and 
we have masslcss states at strong coupling. To do so, wc evaluate the stability of the 
saddle point (3.14). The effective potential for ip is obtained by integrating out a from 
(3.12), 

^eff(<^) = ^(e^'^^'-4V^)- (3.16) 

This shows that we have a second order transition at v — Vc — A."^ , between a broken 
symmetry phase v > and a phase v < where the Z2(iv) is restored. This is the 
"strong coupling" region of parameter space, where the saddle point is 

{A^ 
* rv (3.17) 

(p* = 0. 



24 



It is interesting to investigate in details what liappens at the point v — A?. In the 1/N 
expansion, the effective couphng constant ^^eff goes to infinity at that point, which is 
directly related to the fact that we have massless degrees of freedom. This phenomenon 
is also known to occur in supersymmetric gauge theories when a magnetic monopole 
becomes massless. Actually we will see in Section 4 and 5 that effects which are non- 
perturbative in 1/N make the effective coupling constant finite at the transition point 
V = K^. The critical theory will be argued to be very similar to an Argyrcs-Douglas 
CFT [7], where both magnetic monopoles (topologically stable solitons) and electrically 
charged particles (created by the elementary fields) can become massless. In our case, 
it is clear that the stable SO(iV — 1) singlet soliton found at the end of Section 2.3.3 
must be massless when v — h?. This is simply due to the fact that the two degenerate 
minima of the potential (3.16) merge at that point, and any solitonic solution to the 
effective action must then become trivial. To find out whether any perturbative state 
could become massless together with the soliton, it is natural to first look at the mass 
of the elementary fields $j. One can calculate straightforwardly the two point functions 
in the leading N ^ oo approximation. In momentum space they are simply 

which shows that the fields ^i, 1 < i < N — 1, create free particles of masses = y^aT- 
The mjS thus take their classical values = ^/v as long as f > A^, and then take the 
constant value rrii = A when v < A"^: at strong coupling, a mass gap is created for these 
particles, as the standard lore suggests. To find massless perturbative states ai v — A"^, 
we thus have to look more carefully at the menu of stable particles described in Section 
2.3. As explained there, though the interactions between the elementary particles are 
of order 1/N, a. bound state whose binding energy is of order N^ can be formed because 
of a dramatic stabilization due to the mixing between the different fiavors. This bound 
state is created by the operator J^^^^ but can equivalently be associated with 



1 



N-l 



l-9^J2^l (3.19) 



i=l 



and thus its mass must appear as a pole in the {(f(f) correlation function. By expanding 

= H — ^) a — a^-\ — ^) (3.20) 

Vn Vn 



we immediately obtain from (3.11) the leading term in the large N expansion, 

ixxKp) = rrz — 2 ^ — -2- (3-21) 
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The function s^^^ is defined in Appendix B. The bound state mass rrib thus satisfies the 
equation 

2 

ml = a^-v- ^* — 2 V' (3-22) 

At strong couphng the solution is simply 

mb = VA2 - V, for v < Al (3.23) 

We see that when v = and wc have the full 0(iV) symmetry, mb = A = rrii, and the 
spectrum consists in a O(A^) vector [37] made up by the A^ — 1 elementary particles and 
the (f field. We also discover that when v = A^ the mass of the bound state vanishes! 
When v > A^, we have to solve 



4f V / Av 

^-1 In— = 7r-2 arctanW^-1- (3.24) 

ml y ml ^ ^ 

When V » A^, this yields 

mb = 2^^ (l - 1 {Ng'f + ±- (Ng'f + 0((Ng')\ 1/N)) , (3.25) 



in agreement with (2.51). In the vicinity oi v — A^, (3.23) and (3.24) predicts 

ml ~ A^ (^1 - j for V ^ A^, v < A^, 



ml ~ 2A^ (^-^ - Ij for ^; ^ A^ v> A^. (3.26) 



We have depicted in Figure 3 the masses of the elementary particles and of the 
bound state as a function of y^, as given by the leading term in the 1/N expansion. 
As we will see in section 3.4, this leading approximation is actually incorrect near 
V — A^, and in particular the formulas (3.26) are wrong. However, the most important 
fact that the soliton and the bound state are massless will remain, and we will see that 
the correct form of (3.26) is simply 

ml oc ml^^ (xA^(l-^^ for v Al (3.27) 
3.3 Correlators and quantum corrected metric 

We continue to restrict ourselves to the symmetric case Vi — --- — vi^^im. this subsec- 
tion. This makes the formulas simpler, but do not affect the physics we want to discuss. 
Our aim is to obtain the large N formulas for the 5" matrix, generalizing (2.48), and 
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Figure 3: The mass of the elementary particles created by the $jS, 1 < i < — 1 (dashed 
line) and of the <^-<^ bound state (solid line), as a function of \/v = = • • • = ^/v^^, in 
the leading 1/N approximation. 



to compute the quantum corrections to the metric on the sphere target space. The 
generating functional Z[J] for the correlation functions is given by a path integral 



Z[J] 



VxD 13 exp 



-(iV - 1) Seff[a, + /3/ViV, + x/ViV] 

N-l 



+ - / d^x ^ Ji 



i=l 



+ a, + 



Ji 



VxV(3exp -{N-l) Sefr[a* + (3/VN, ip, + x/VN] 



(3.28) 



The large Feynman diagrams contributing to the four point function (^jj^jj^ia^M) 
are indicated in Figure 4. The wavy lines correspond to the (5 field propagator 

{mip)= .,^ ^'^''^J^ , 2- (3-29) 

The functions Si, S2 and S3 defined by (2.47) can then be deduced from the four- 
point function. At weak coupling {v > A^), we have 

'^1 = FT^= 1 I /\r 2 — ^75171 ^' '^2 = 1, S3 = Si{s^u), (3.30) 
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Figure 4: Large Feynman diagrams for the four-point function. The wavy hne corresponds 
to the f3 field propagator, which is itself a sum of bubble diagrams in ordinary perturbation 
theory. 

or equivalently in terms of the rapidity > such that s = 4f cosh^(0/2) and u = 
-iv smh^{9/2), 

^iW = ^^-r^ ^ir-JT ' ^2^ = 1, Ss{9) = S,{z7r-9). 



8.^ tanh(^/2),^^^.^ 



ITT 



^'^^47rtanh(0/2) 

(3.31) 

is defined as usual to be the low energy coupling, i.e. the running coupling (2.27) 
evaluated at = ^/v. These formulas illustrate nicely how the expansion works 
and produces non-perturbative results. The genuinely non-perturbative information 
lies on the full 9 dependence of the denominators, not on the fact that appears 
in these denominators. In particular, when 9 ^ oo, the non-perturbative corrections 
become dominant, even if f ^ A^, but due to asymptotic freedom, we know that in 
this regime the S matrix can also be deduced from perturbation theory by using the 
renormalization group. Indeed we have 

S^{9) = ^ = (3.32) 

which is the perturbative result (2.48) at high energy, but with the coupling g evaluated 
at the center of mass energy, which is different from the low energy coupling. In the 
strong coupling region f < A^, we have 

i 1 zA^ 

■^l = g 2 AT I ~^V£U~2 X^' "52 = 1 , ' 53 = 5i(s^u), (3.33) 

or equivalently 
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The functions SiS are independent of v in this regime, but this is true only in the 
leading approximation. The formulas (3.33) show that the $jS, l<i<A^ — 1, no 
longer form a two particles bound state for v < A^, and thus that the field must be 
considered as an independent degree of freedom in this region. Clearly the geometric 
interpretation of the $iS as living on a sphere cannot be valid for v < A^. To better 
understand what is really going on, it is natural to compute the quantum corrections to 
the metric on the target space sphere as a function of v. This amounts to computing the 
coefficient of in a low energy expansion of all the connected, one-particle irreducible, 
2n-point functions ■ ■ • $i2„)- The SO(iV— 1) symmetry restricts a priori the metric 
to be of the form 

N-1 N-1 

= /i (E ^i) + /2 (E ^i) ^^^^ (3-35) 

k=l k=l 

for unknown functions /i and /2. Actually, we are going to show that in the leading 
1/N approximation and for v > 

gij-hi^ ' (3.36) 

k=l 

which is simply the SO(A'") invariant metric on a sphere of radius 1/geS- This very 
simple result is equivalent to the statement that the terms in the low energy ex- 
pansion of the 2n point functions is simply given by the perturbative, tree level result, 
with the bare coupling being replaced by the low energy coupling Qcs- This simpli- 
fication occurs because the low energy expansion of the j3 field propagator (3.29) is 
{(3(3){p) = —Ngl^p^ + O(p^), and thus only one such propagator can appear in a dia- 
gram contributing to the metric. A connected diagram of this type for the 2n points 
function must then be entirely constructed from one {iSjd) propagator, n — 2 vertices 
l3-)(^/\fN coming from the expansion of Seflf (3.11,3.20), and 2n — A contractions between 
/3 and x which each give a factor of -s/Ngl^, in addition to the n ^^P/y/N vertices. 
Such diagrams at low energies depend on N and g^s only through a multiplicative fac- 
tor Ar^2^ Ar-("-2)/2 (Ar^2^)(2"-4)/2 7v-"/2 = N^-r^ (^V^lff)""^ = 9es'^^- This dependence 
on the coupling is the same as the one found in a tree level perturbative calculation. 
The dependence in the momenta must then also coincide, since the large N result must 
reduce to the tree level result when g^s <^ 1. This proves (3.36). 

We thus see that the effective target space sphere shrinks symmetrically to a point 
when i> — > A^, i> > A^. The geometrical interpretation of the a model is thus lost in the 
strong coupling region v < A^, which corresponds in some sense to a continuation of the 
sphere to negative radius. Similar continuations have already been seen in the context 
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Figure 5: Typical diagrams relevant to the calculation of the metric and contributing to the 
6-point (left) and 8-point (center and right) functions. The {(5(5) propagator is represented 
by a wavy line, and the (/?x) contraction by a dashed line with an arrow pointing toward the 
X insertion. 

of supersymmetric models [39], and when the a model is conformal they actually play 
an important role in understanding the topology changing transitions in string theory 
(see for example [40]). We will further discuss the consequences of this interpretation 
in Section 5. 




Figure 6: The effective target space sphere shrinks symmetrically to a point when v — 
in the leading 1/7V approximation. At strong coupling, the geometrical interpretation of 
the cr model is lost. This is reminiscent of the stringy geometry described in the context of 
supersymmetric non-linear a model. 

3.4 The critical hypersurface 

Let us now discuss the general case where the Vi>Q can be distinct. Formulas (3.11) — 
(3.14) and (3.17) are replaced respectively by 
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a-h ^ N-1 a-h- 




87riV " eh? 

1=1 



^ ^ a^ — hi (3.39) 



A/" ^ 



TV ^ A2 

i=l 

Weak coupling : { 2 _ 1 ""i (3.40) 

E — Al; _ 
^..^o . ^ A2 ~ (3.41) 

= 0. 

Note that the various sums over i should be more rigorously replaced by integrals with 
the help of (3.2). This will always be understood in the following. The border between 
the weak coupling region ($Ar) 7^ and the strong coupling region {^n) = is a 
codimension one hypersurface one which both the soliton and the bound state become 
massless. The equation for this quantum hypersurface of singularity 7iq in Aljv is 
immediately deduced from (3.40) and (3.41), 

N-1 

Hc^HMn ■■ l[vi^A''^''-'\ (3.42) 

i=l 

The critical point vi — • ■ ■ — vn-i = Vc = discussed previously is of course on TYq. 
Tiq has also other sheets in the other, physically equivalent, regions Alj of parameter 
space (2.43), whose equations are simply 

n^nMi : = A2(^-^). (3.43) 

We thus get a preliminary picture (it is only partially correct as will be discussed in 
Section 4) of the quantum space of parameters, in the large limit. 

It is very likely that the qualitative picture is still valid even for values of as 
small as A^ = 3 (we will give strong evidence that this is the case in Section 5), and we 
have used equations (3.42) and (3.43) in the cases N — 3 and A^ = 4 to make Figure 7. 
This figure should be compared with the classical case given in Figure 1. One might 
have expected that the singularities on parameter space would have disappeared in 
the quantum case, a mass gap being created, and the 'Z2{n) symmetry being restored 



31 




Figure 7: A preliminary picture of the quantum space of parameters A^q in the cases N = 3 
(left) and A" = 4 (right). 

through a first order phase transition. We see that on the contrary the classical singular 
hyperplanes split into two in the quantum case, and the 1'2{n) symmetry is restored 
through a second order phase transition. We discuss in the next subsection how we can 
understand the nature (Ising) of this transition in the large A^ limit, and we will also 
discuss it at finite A^ in Section 5. 

We can also write down the formulas generalizing (3.23) and (3.24) to the case of 
arbitrary positive ViS. At strong coupling, the mass mb satisfies 

= Va*- hN, (3.44) 

with a* determined by (3.41), and at weak coupling we have 

N-l N-l ^ ^ 

V In ^ = 2 V , arctan , • (3.45) 
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3.5 Beyond the 1/N expansion: the Ginzburg-Landau description 

In the vicinity of the critical hypersurface Tiq, the usual expansion is inconsistent 
because of infrared divergences. These divergences are due to two different effects. 
First, the % field becomes massless and thus the (xx) propagator (3.21) goes to 
Consequently the 1/N expansion cannot be defined because of the usual IR divergences 
associated with massless scalar propagators in two dimensions. Second, and indepen- 
dently of the fact that we are in two dimensions, the 1/N corrections are calculated by 
taking into account interaction terms corresponding to relevant operators. Such inter- 
actions always produce IR divergences in perturbing around a massless theory. This 
is strictly analogous to the textbook discussion of the corrections to the mean field 
approximation for the (f)^ theory near dimension 4, except that we are dealing presently 
with the perturbation theory instead of the ordinary perturbation theory. To 
analyse the divergences, one could study in details the large N Feynman graphs and 
isolate the most divergent contributions. Instead, we will proceed in a slightly more 
heuristic, but actually equivalent, way. In a first step, we would like to integrate out 
the field a = a* + P/\/N, since the divergences are due to the masslessness of x- At 
large N, this can be done by solving the variational equation 



which at low energy reduces to 



dVeS 



By using (3.38) and (3.37), we then obtain a low energy effective action valid in the 
vicinity of Tiq, and depending on x = only, 

^eff [X] = / d'x d^xdaX - Y + ^ + Oix'/N')) , (3.48) 



where we have defined 



and 



1 1 

- = ^— y -, (3.49) 
V N-1 ^ V, ^ ' 



N-1 
i=l 

The equation for Hq is simply Sv — 0. Higher derivative corrections to (3.47) produce 
terms like X^daX^aX/^ which are irrelevant in the IR. Similarly, we have dropped in 
(3.48) irrelevant terms or order x^^ /N^~^, k > 3. In a 1/N expansion at Sv — 0, their 
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contribution is subleading with respect to the x^/^ term, to a^/ orders in 1/N. Indeed, 
IR divergences compensate for the factors coming from the x^/^ interaction, 
which should thus be considered as giving leading contributions when Sv ^ 0. A 
way to understand this is to rescale the space-time and momentum variables 

x^^/Nx', p = p'/^/N, (3.51) 

which eliminates any N dependence in (3.48) on Tiq, 

SM = / rfV (^^dUdU + ^Vx' + OixVN)). (3.52) 

The IR properties p — > can now be studied by working at fixed p' — py/N and 
using a simple minded 1 /N expansion in the p' variables. This amounts to perturbing 
around the interacting theory (3.52) where the term is taken into account exactly, 
and contrary to the perturbation around the free massless theory this is perfectly well- 
defined. 

We recognize in (3.48) the Ginzburg-Landau description of an Ising critical point 
[41]. The low energy physics on TYq {6v = 0) is thus governed by an Ising CFT. Formula 
(3.27) is then simply derived by noting that a deviation from the critical value v — 
corresponds to turning on the energy operator in the Ising CFT, a dimension one op- 
erator. Note that displacements along the critical surface correspond to operators that 
are irrelevant in the IR, since the Ising CFT docs not have any marginal deformation. 
The appearance of a non-trivial interacting CFT at a point where both a soliton and a 
pcrturbativc state become massless is of course very reminiscent of an Argyres-Douglas 
CFT [7], though in our case it appears to correspond to a very banal critical theory 
with a local description (3.48). Nevertheless, we will see in Section 5 that a natural 
way to understand the appearance of this interacting CFT, independently of a large N 
approximation, is through a description in terms of a more exotic non-local theory of 
electric and magnetic charges (spin waves and vortices) . The Ising CFT is then nothing 
but a strong coupling fixed point of this non-local theory. This latter description is 
more akin to what one may expect in the context of gauge theories, where simple local 
description a la Ginzburg-Landau may not exist. 

The reader may be concerned at this point by the fact that even though the stan- 
dard 1/N expansion is not valid in the vicinity of 7iq, the very existence of TCq itself 
has been demonstrated in this simple framework. To investigate the consistency of the 
analysis, one should compute the 1/N corrections to the equation for Hq (3.42), and 
see if they are small. This calculation, which illustrates several properties of the 1/N 
expansion, is presented in Appendix C. The result is that though we do encounter IR 
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divergences when computing the corrections to (3.42), they are mild logarithmic diver- 
gences. The only consequence is that the leading corrections, instead of being of order 
1/N, are of order (In N)/N. This is still a small correction when — > oo, as required. 



4. The quantum space of parameters Ad^ 

The equation (3.42) for Hq is not reliable when a small number of the ViS goes to zero. 
Indeed, this equation is obtained in the leading N ^ oo approximation in the form 

N-l 

n,nM^ : ^E^^x^^o- (4-1) 

i=l 

Then, when for example vn-i is much smaller than the other ViS, the corresponding 
term in the sum (4.1), which is of order 1/N, inconsistently becomes larger than the sum 
of the other terms, which is of order N^. Naively, one might expect that this problem 
could be eliminated by taking into account the corrections to the equation for TYq 
(C.30), but this is not true: the fact that the large N counting is changed in the limit 
vtv-i — > affects all orders of the 1/N expansion used in Section 3. 

4.1 Weak coupling 

Actually, we can readily understand that (4.1) is not correct, even qualitatively, on the 
whole parameter space. Indeed, when fjv-i — and all the other ViS are much larger 
than A^, the low energy physics is governed by a weakly coupled 0(2) non-linear a 
model. This is a massless phase, and thus we see that at weak coupling the classical 
(Figure 1) and quantum hypersurface of singularities must coincide. This is unlike the 
prediction based on (4.1), see Figure 7. More precisely, for v^q^i <^ Vi and Vi ^ A^, 
1 < i < — 2, the low energy effective action is reliably determined by a one-loop 
perturbative calculation (2.24) to be 

^eff = ^ J d^x (djdj - ^ COS29') (4.2) 

with 

gls 47r A2(^-2) 

As already discussed in Section 2, we get a sine-Gordon model whose coupling in 
standard normalizations is 

PsG = 2^efr. (4.4) 
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It is important to note that, though TYq = 7ic\ at weak couphng, the classical and 
quantum low energy physics are drastically different. Classically, we have a free massless 
theory independent of Qes- Quantum mechanically, we have the non-trivial CFT of a 
boson compactified on a circle of radius 1/ g^s- Even at very small coupling, the physics 
(anomalous dimensions, etc. . . ) depends on ^^eflf in a non-trivial way. The fact that 
the classical and quantum singularity structure can coincide though the corresponding 
physics arc different is also seen in supersymmetric gauge theories, for example in A'' = 1 
super Yang-Mills with A^f = A^c + 1 [2]. 

The description in terms of (4.2) is very useful at small g^s (for a review on the 
sine-Gordon model including the original references, see e.g. [42]). We can for example 
use (4.2) to compute the mass of the two particle bound state, 

A^^i ( gl^^ gl^ 

— 



& V 2nJ 2(1 -gy 2^) 

This formula should be compared to the analogous formulas (2.51) and (3.45) which 
were obtained in different regimes. The mass of the soliton (2.55) is also known, 

m.. = ^fl-^). (4.6) 



The action (4.2) also predicts that the bound state becomes unstable when g'^g > 27r/3, 
and that we have a massless phase for g'^g > 2ti even when Vi ^ 0. This is of course 
incorrect, and (4.2) is not a good description of the physics at strong coupling. We 
know that the non-abelian degrees of freedom must come into the game and create 
a mass gap A. Moreover, the results of Section 3 indicates that for vn-i 0, there 
should be a critical value of g^s at which both the bound state and the soliton become 
massless and the low energy theory is an Ising CFT. This is consistent with the idea that 
(4.5) only gives an upper bound on the mass of the bound state, which can actually 
be significantly lower when the coupling increases due to the mixing with the other 
degrees of freedom (see Section 2). 

4.2 Strong coupling 

To investigate the strongly coupled region, the idea is to modify the large N expansion 
of Section 3 in order to treat exactly the fields and ^n-i is the region of small 
I'jv-i- This amounts to integrating out the N — 2 fields . . . from (2.31), while 
keeping explicitly $jv = y/N^p and ^n-i — y/Nip'- The effective action, which replaces 
(3.37), is then simply 

Seff[Qi, ^,^\ = j d X\- da^pda^P + - O^P + P H (f - 1 
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i=l 

We can then repeat straightforwardly the arguments of Sections 3.4 and 3.5. The large 
N low energy effective action at small vn-i coincide with the one-loop result (4.2) 
except in a vanishingly small region in parameter space corresponding to 5v <S 
where 

i=l 

with 

1 1 1 

't=l 

In this region, the fluctuations of the radius of the effective target space circle become 
important, and the low energy physics is governed by the Ginzburg- Landau action 

^efr = / d^x (i d^xdaX + \ + ^ X'^ - y + x'^) + ^ (x^ + x'^)') • 

(4.10) 

This type of Ginzburg-Landau theory has already been studied [43], following the ideas 
of [41]. If V]ss-\ = and Sv > 0, it gives a description of the compactified boson CFT. 
The marginal operator corresponding to a change in the radius is x^ + x'^- We thus get 
a smooth interpolation with the physics described by (4.2). If v^^i > 0, we have an 
Ising critical point corresponding to the restoration of the 1'2{n) symmetry at 5v — 0, 
which means that the equation for Hq is 

iV-2 

H^HMn ■■ Y[vi^A^^''-^\ for < < A'. (4.11) 

i=l 

We thus also obtain a smooth interpolation of the physics described by (3.48), but we 
now see that Hq must intersect the hyperplane fAr_i = 0. Similarly, if Vn^i < 0, we 
have an Ising critical point corresponding to the restoration of the 'I^2{n-i) symmetry 
at 5v — vn-i, and thus 

N-2 

HqnMN-i : Ylvi= (A'e"^-i/^) , for - A' < vn-i < 0. (4.12) 
1=1 

For vjsr^i — and Sv — 0, the two Ising CFT found for i'at-i > and i'at-i < 
respectively are coupled through an 0(2) invariant interaction, giving an Ashkin- 
Teller critical point. This CFT is equivalent to the compactified boson CFT described 
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by (4.2) for vn-i = at g'^g = 7r/2, which is itself equivalent to the Z2 orbifold 
CFT at the self-dual point [44] (for an elementary discussion, see [45]). The marginal 
operator corresponding to the variation of the orbifold radius is x^x'^ [43] , but this is 
not associated to any microscopic operator in the Neumann model. For this reason, the 
equivalence with an orbifold theory does not seem to play a particular role in our model. 
What is extremely significant, however, is that the value g'^f^ = 7i/2 also corresponds to 
the critical coupling for a Kosterlitz-Thouless phase transition [46]. For f7v_i = and 
6v < 0, (4.10) is indeed believed [43] to describe a massive phase, and no longer the 
compactified boson CFT. This is a very interesting physics: the creation of the mass 
gap in our model, which is fundamentally due to the non-abelian nature of the degrees 
of freedom, can be understood as coming from a condensation of vortices in the abelian 
description (4.2). This aspect will be discussed in great details in the next Section. 

The fact that the transition must occur at g^f^ = 7r/2, and not at g^Q = 00 as 
predicted by the N ^ 00 approximation (4.11, 4.3), is an exact result that cannot be 
deduced in any simple approximation scheme. Physically, it means that the effective 
target space does not literally shrinks to a point, but that the effective radius never- 
theless becomes so small that its quantum fluctuations are important and any classical 
geometric interpretation of the target space is impossible. This is likely to be true any- 
where on Hq, and thus for example the endpoint of Figure 6 should rather be viewed 
as a very small, but non-vanishing, fluctuating quantum sphere. 

4.3 The ansatz for the equation of Hq 

We would like to present an equation for TYq valid in the large N limit on the whole 
parameter space Ai, that interpolates smoothly between the different regimes (3.42, 
3.43, 4.11, 4.12). To get such an equation, we clearly need to treat all the regions Aii 
(2.43) symmetrically. It does not seem to be possible to do this rigorously, for reasons 
explained at lenght above: we need different 1/N expansions in the different regions Mi, 
and yet another expansion to understand the transition between the different regions. 
However, there is a fleld, the Lagrange multiplier a, that we can treat symmetrically 
on the whole parameter space. In particular, we can compute an A'^ — > 00 effective 
potential for a by integrating out the flelds $i from (2.31), 




(4.13) 



i=l 



The saddle point a* satisfles 




N-l 



(4.14) 



i=l 
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Figure 8: The quantum space of parameters Aiq in ttie cases N = 3 (left) and N = A (right). 
We have chosen arbitrarily A/ A = 3/2. 

It is easy to see that (4.14) is consistent with our preceding formulas (3.40, 3.41). For 
example, if we consider the region A^jv at weak coupling (A small), (4.14) implies an 
expansion a* — + ^^^^'^^11^=1^ Vi + • ■ •, and thus the corrections to (3.40) go to 
zero except on a vanishingly small region of AIat when — > oo. The advantage of 
(4.14) is that it is smooth over the entire A4. An improved equation for TYq in AIat is 
then to use (3.39), 

N-l 

n(«*-M=A2M, (4.15) 
1=1 

but with a* now being determined by (4.14). A natural generalization of (4.15) on the 
whole of parameter space is 

N N 

i=i j=i 

where we have introduced a "phenomenological" constant k > \ which, though it can- 
not be seen in the N ^ oo limit, is necessary to make Tiq intersect with the hyperplanes 
separating the different regions A4i. k can be considered as a non-perturbative quantity 
analogous to the Kosterhtz-Thouless critical coupling g^s — 7r/2 that cannot be calcu- 
lated in a large A^ expansion. The best justification we can find for (4.16) is by looking 
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a posteriori to its consequences: by construction, is reproduces the correct result in the 
regions Mi separately (Section 3.4), and it can be checked easily that it also gives a 
correct interpolation between the different regions (Section 4.2). Though our results 
have been obtained at large A^, it is very likely that they are qualitatively valid even 
for N — 3, and we have used (4.14) and (4.16) to obtain a sketch of the quantum space 
of parameters Aiq in the cases = 3 and = 4 in Figure 8. 

In our model, it is necessary to cross a surface of singularities to go from weak 
coupling to strong coupling. This is of course not true in general. For example, one 
may turn on a magnetic field (a {v, 2) mass term in the terminology of Section 2) in 
addition to the (s, 2) term considered so far. This would make the phase transitions 
first order. 

5. The non-local theory of electric and magnetic charges 
5.1 General discussion 

The picture of the quantum space of parameters in Figure 8 clearly shows that we 
have two quahtatively difTerent regimes in our model. Outside the hypersurface of 
singularities Hq, the usual geometrical interpretation of the non-linear a model as a 
sum over maps from to S^"^ is valid. The use of the coordinates on the sphere 
as the fundamental physical degrees of freedom is then justified. On the contrary, 
in the interior of Tiq, the geometrical interpretation breaks down. It is of course still 
possible to define formally the model as a sum over maps from to S^^^, but the path 
integral will be dominated by highly singular maps from which the image of a smooth 
target space manifold cannot be reconstructed. A very interesting possibility is that 
there could be a (non-local) change of variables in the path integral to new degrees 
of freedom that would give a natural description of the physics at strong couphng. 
These new degrees of freedom may, or, more probably in our case, may not, have 
a geometrical interpretation in terms of coordinates on a new target space. In the 
case oi H — 2 supersymmetric non-linear a models, this change of variables can be 
argued to exist in many cases and is called the "mirror map" (see e.g. [40], [47]). In 
the case of the integrable O(A^) non-linear a model, which corresponds to the origin 
of our space of parameters (Figure 8), it is actually a long standing problem to find 
these "good" degrees of freedom (see e.g. [48]). When A^ = 3, a possibility is to use 
variables associated with instantons. In our model, instanton calculations can be done 
reliably at weak coupling since the ViS provide there an IR cutoff. However, I think 
that it is unlikely that instantons, which are smooth configurations of the original fields, 
can account for the physics near and in the interior of Tiq (it would be interesting to 



40 



investigate this point, though). Moreover, one would hke to have a unified description 
for all > 3, and not only for = 3, since the physics is qualitatively the same for 
all those cases. It is more natural to suspect that the relevant degrees of freedom at 
strong coupling are more singular configurations. Figure 6 suggests that configurations 
of the fields where 

$ = (5.1) 

at some space-time points will be conspicuous. Around such points, the angular coordi- 
nates on the sphere can vary abruptly. These configurations are suppressed classically 
by the kinetic energy term in the action, but at strong coupling we see that they can 
actually contribute significantly to the path integral. For example, in the case of the 
two-sphere, we have two coordinates 6 and with the identifications 

(^,(/)) = (^ + 27r,(/)), (e,0) = (^,(/) + 27r), = + tt), (5.2) 

which suggest that relevant degrees of freedom could be vortices in the 9 and (j) vari- 
ables as well as configurations where 9 has square root branch cuts, as in a Z2 orbifold. 
Unfortunately, I do not know how to take into account this three types of defects to- 
gether in an 0(3) invariant way. Nevertheless, we can more modestly look at the regime 
vm-i — > studied in Section 4 where the low energy effective theory is abelian. When 
vn-i — 0, we have demonstrated at large N that the transition from weak (massless 
phase) to strong (massive phase) coupling is of the Kosterlitz-Thouless type. It is of 
course well-known [46] that such a transition can be understood as being triggered by 
a condensation of vortices 9 = 9 + 27r. Let us stress an important non-trivial point of 
principle here. The low energy action at vn-\ = 0, 

-^eff^^ [d^xdjdj, (5.3) 

does not by itself predicts a Kosterlitz-Thouless transition at g^f^ = 7r/2, but rather a 
massless phase at any g^s- The action (5.3) is often used as a continuum version of 
the 0(2) spin model defined on a lattice, g^g playing the role of a temperature. The 
lattice model is the object of interest in condensed matter physics, and was studied 
by Kosterlitz and Thouless [46]. In the lattice hamiltonian, configurations for which 
9 9 + 2-11 abruptly are not suppressed; the suppression is an artefact of the continuum 
description given by (5.3). By introducing the vortices by hand in the continuum 
formulation (which means that we choose a non-zero fugacity for the singular field 
configurations corresponding to vortices in the definition of the path integral) , one can 
get a good modehng of the lattice model [46] , and in particular predict the mass gap at 
high temperature and the correct behaviour at the transition. In our case, however, we 
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are not dealing with lattice models, but with a bona fide field theory. What has been 
proven in the large N limit in Section 4 is that a non-zero fugacity for the vortices is 
dynamically generated in the model because of the non-abelian (continuous) degrees of 
freedom. The physical origin of the condensation of vortices is thus completely different 
in our case and in the well-known abelian lattice model. Note that the vortices are 
not instantons of the non-abelian model and can exist for all N. This suggests that, 
more generally, the O(A^) model may be naturally described in terms of "non-abelian" 
generalized vortices (5.2). 

In the following, we are going to study how the standard picture of vortex con- 
densation is modified when we go to the regime where wat-i is small but non-zero. In 
Section 3 and 4 we obtained in the large N limit a simple Ginzburg-Landau description 
of the physics on Tiq. We do not expect to get as simple a description in the present 
context. In general, both weak coupling and strong coupling degrees of freedom can 
be relevant in the vicinity of Tiq, where the transition between the weak coupling and 
strong coupling behaviours takes place, giving an unusual non-local description of the 
low energy physics. The same kind on physics has of course already been studied in 
condensed matter works on abelian lattice models. A cos p6 potential in the lagrangian 
is interpreted in this context as a p-fold symmetry breaking term modeling crystal 
anisotropy. When p > 4, a perturbative analysis a la Kosterlitz is possible ([49], [50]), 
but when p < 3, and in particular in the case p = 2 we are interested in (4.2), this is 
not possible as we will review. Even in those cases, the nature of the physics is rather 
clear from the lattice point of view in particular (see the Appendix of [51] or [52]). The 
general ideas presented later in this Section are thus certainly not new, though a com- 
prehensive analysis of the kind we offer does not seem to have appeared. Our main goal 
is really to emphasize the striking similarity between our non-supersymmetric system 
and the physics in the vicinity of an Argyres-Douglas point in supersymmetric gauge 
theories. Whether similar things can happen in non-supersymmetric four dimensional 
theories is of course still an open problem. 

Before we turn to these points, let us make a simple heuristic remark. In our model, 
we have a clear interpretation of the fact that a compactified boson can be equivalent 
to the 0(2) symmetric Ashkin- Teller model, and that this happens precisely at the 
Kosterlitz-Thouless coupling. One may then wonder whether a natural explanation 
for the fact that these CFTs are also equivalent to the Z2 orbifold at the self-dual 
radius could emerge. A naive argument in our context could be the following: since the 
Kosterlitz-Thouless transition is triggered by the non-abelian degrees of freedom, the 
angle 9 should know somehow that is lives on a sphere, and an or bif old-like identification 
9 = —9 should then be implemented (5.2). However, the orbifold of the circle theory 
at the Kosterhtz-Thouless radius is not the orbifold at the self-dual radius. The naive 
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interpretation thus does no make sense as it stands. It would be interesting to try to 
include the effect of the other $ coordinates on the sphere to see if this picture might 
be made consistent. 



5.2 The non-local description in the abelian regime 

A general configuration of the field 6, including vortices, can be written 

e{x) = eo{x) + e,{x) (5.4) 

where is the non-singular part and 

Ov{x) = '^niaig{x - Xi) (5.5) 

i 

corresponds to vortices ("magnetic charges") of charge rii centered at the points Xi. 
When the total charge ^ rii is zero, the action (4.2) can be finite and is given by 



I Ob Ob j I 



1 /" TT . 

Ses{ni]Xi) = / (PxdaOodado o- niUj In 

^/^=4n(|^)^--n(f^) 



.(5.6) 



We recognize the standard Coulomb interaction between magnetic charges in two di- 
mensions, but we have also a complicated interaction term between the magnetic 
charges and the fundamental spin waves ("electric charges"). We have noted Zi = 
x} + ixf and Zi = xj — ixf. The charges rii = ±1 are the most relevant, as can be 
checked straightforwardly along the lines of the calculation presented below. We will 
thus write the partition function in term of the dynamically generated fugacity F = /A^ 
(/ is a dimensionless constant) for the rii — ±1 vortices as 

z^J2 -7:^ /n(^'^A) vOoe-'^-^'^-'^-y^K (5.7) 

n=l ^ '■' i=l 

where the XiS and yiS parametrize the positions of the charge -|-1 and charge —1 vortices 
respectively, and Ses[Oo; Xi, i/i] is the action (5.6) for such a configuration. By expanding 
the exponential in powers of i'at-i and calculating the path integral over 9q, we get 
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^= E ' " faff' 'faii /n(<<vH)n(''=-.A<) 



J]^ {\Up - Uq\\Wp - Wqlf^"^^"" Y\ \Up - Wq^ 



l<i<jr'<n 

TT (^p - (Wp - Vt) jUp - 7/^) (Wp - Xj) 
l<p<m 

The second line in (5.8) corresponds to the standard mass expansion of the original 
sine-Gordon model (4.2), the third line corresponds to the interaction between vor- 
tices, which is also described by a sine-Gordon model of the type (4.2) but with a 
coupling gigg = ir/geff, and the fourth line corresponds to the interactions between the 
electric and magnetic charges. This partition function is reproduced, including the 
electric/magnetic interactions, by the perturbative formula 



-Sint\ 

/ free 



(5.9) 



where the average is computed with the free field weight (5.3) and 

Si^t = - fd'^xi cos 26 + 2/A2 cos ^ 6] . (5.10) 

J V4^efr ^eff J 

The dual field 9 is defined by 

d9 = -i*de. (5.11) 

The formula (5.9) for Z is perturbative, since it involves an average over free fields. 
Interestingly, an off-shell, non-perturbative, formulation of this type of theories exists 
[53]. The basic idea is to treat the fields 9 and 9 as independent variables, and try to 
find an action whose equations of motions give (5.11) when the interactions are turned 
off. By writing (5.11) in the minkowskian, 

do9 = -di9, di9 = -do9, (5.12) 

we see that di9 is in some sense the canonical momentum associated with 9. The phase 
space path integral then suggests to try the action [53] 



d'^x (do9dJ + do9d^9 + {d^9f + {di9f^ + 5int[^, ~9] (5.13) 
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which indeed yields the correct equations of motion. If S'int depends on only, then by 
integrating out 6 we would recover the standard Lorentz invariant action for a scalar, 
but in general it is impossible to write down a manifestly Lorentz invariant action 
for both 6 and its dual. Operators associated with a dyon of electric charge and 
magnetic charge are 

We thus see that our model (5.10) corresponds to the case where both a magnetic 
charge (the vortex) (ng = 0, rim — 1) and an electric charge (rig = 2, rim — 0) are 
included. By rescaling the fields 

9^gef£(f), 9^gef£4> (5.15) 

and defining 

Fm = 4fA' (5.16) 

we can bring the action in a form 

Sera = / {do<pdi4> + do4>d,<P + {di<pf + {d,4>)^ + F, cos{2g,s<p) + Fm cos(27r0/(?efr)) 

(5.17) 

which is manifestly invariant under a strong/ weak coupling S duality, 

fl'eff^^ — ' Fe^^Fm, 0. (5.18) 

Note that S is not the usual T duality that would exchange the unit magnetic charge 
with the unit electric charge and g^s with 2-11 /gcs_^ in the same way as the monodromy 
at an Argyrcs-Douglas point is not the full duality group. 

The quantization of (5.17) can be done in the standard way. We have two second 
class constraints, 

n = -lai(^, n = -lai(/., (5.19) 

and the Dirac bracket yields the equal time commutation relation 

^iQ{x^-x'^) (5.20) 

where is the Heavyside step function. We see that the non-locality of our theory, 
which describes the interaction of electric and magnetic degrees of freedom, is simply 
encoded in this commutation relation. Let us note finally that actions describing both 
electric and magnetic charges in four dimensions, and that should play a role in un- 
derstanding Argyres- Douglas CFTs, have been constructed along the lines of (5.17) 
[54]. 
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5.3 The strongly coupled fixed point a la Argyres-Douglas 

The standard way to treat (5.17) would be to use perturbation theory in and 
([49], [55]). This makes sense when the corresponding operators are nearly marginal. 
What are we expecting in our case? Prom the results of Section 4, or from the intuition 
gained in lattice model, we would like to show that (5.17) has an Ising fixed point. If 
we have only one fixed point, which is likely, it must occur at the ^'-invariant point 

S-eff = 5* = VtT, Fe = Fm. (5.21) 

On the other hand, the perturbative dimensions of the "electric" and "magnetic" op- 
erators are respectively 

^elec — ' ^mae; — ~9~' (5.22) 

which would correspond at the fixed point to 

^elec* ~ ^mag* ~ 1) (5.23) 

values far below the threshold of marginality A = 2 where perturbation theory is valid. 
This means that the fixed point we are looking for must be at strong coupling. This 
is also what is expected in four dimensions, and in general it would imply, even in 
two dimensions, that the theory in not tractable. Luckily, in the very particular case 
at hand, we have a powerful tool at our disposal that can solve the problem: we can 
fermionize and replace the scalar by a Dirac spinor ip. Fermionizing is of course a 
very natural thing to do since we are looking for an Ising critical point. Moreover, the 
self-dual dimensions (5.23) strongly suggest that the electric and magnetic operators 
can be viewed as fermion mass terms near the fixed point. Fermionization is a bit 
unusual in our case, however, because the "topological" symmetry — > -|- constant is 
broken by the interaction term in (5.17) and thus we do not have a conserved fermion 
current j oc *d(l) (in other words, the relations (5.11) or (5.12) are no longer correct 
due to the interactions between and 0). The perturbative reasoning of Coleman 
[56] can nevertheless be reproduced straightforwardly, starting from (5.9). In the same 
way as the sine-Gordon potential introduces Dirac mass terms, the cos{2n(j)/ges) term 
introduces Majorana mass terms. More precisely, we can identify 

II COs(25reff0) = -TTIpIp, jl COs{2tT^/ Q^s) = -TT'ljj'^lp, (5.24) 

where we have introduced an arbitrary renormalization scale /i. The theory (5.17) is 
then equivalent to the fermionic theory with lagrangian 

-5ferm = J d?X (l^P-f^d^^P + ill + 1^ ^ _ G ^^^^ ^^^^^ ^^ ^5) 
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with 



(6.26) 



It is useful to decompose the Dirac fermion in terms of two Majorana fermions 
A = — A*" and x = — X*^ such that 



^- \ A + ix 1 / A_ + ix 



In terms of these new variables, and by introducing the light-cone coordinates — 
± and 

M="'^'-^-). M' = ^^<^l±iy. (5.28) 

we have 

'S'ferm = J d^x (^A+9_A+ + x+d-X+ " A_9+A_ - x-d+X- 

-iM' x-K - x+ A_ - "^^^^f^ A_x- A+X+) ■ (5-29) 



We thus see that at the self-dual point (5.21), for which M = 0, we have at low energy 
a free Majorana fermion (A_, x+), which indeed describes an Ising critical point as was 
to be shown. The 5" duality (5.18) acts as 

mJ^-M, M'^M', x±^TX±, A±^A±, (5.30) 

and thus reduces at low energy to the Kramers- Wannier duality of the Ising model. 

I would like to conclude this Section by giving a more rigorous derivation of (5.25) 
and explaining at the same time a point that might puzzle the reader. The puzzle is 
the following. The duality (5.18) should be valid not only at the fixed point, but also 
for any coupling ges 7^ g*- Using (5.18) and (5.30), the four-fermions interaction term 
in (5.25) naively transforms under S as 

S G - - 
Gipj''ijijj^iji — » , ^/ ^7^^^7/^^- (5-31) 

i + G/TT 

This is clearly inconsistent, since by repeating the transformation (5.31), we could make 
the coupling arbitrarily small. The subtlety comes from the definition of the fermion 
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current from which the four-fermions interaction is constructed [57]. The canonical 
definition, that must be used for example in a Hamiltonian approach, uses a point- 
splitting regularization at equal time. However, such a definition does not produce a 
Lorentz vector. A correct definition of the current must actually involve both a space- 
like and a time-like splitting. There is a one-parameter family of Lorentz covariant 
current that can be defined in that way as we will review below. Products of fields taken 
at time-like intervals depend on the dynamics of the theory. The Lorentz invariant four- 
fermions interaction must then depend in some hidden way on the coupling constant 
G, and this invalidates (5.31). Since we don't want to repeat the tedious arguments 
of [57], we are going to show much more straightforwardly how this comes about by 
giving at the same time a rigorous non-perturbative derivation of (5.25) starting from 
(5.17). We will consider the hamiltonian in the Schrodinger picture, 

2 / ,7\ 2 

(5.32) 



\~r] + {'T~] + COs(25feif0) + Fm cos(27r0/5feff) 



1 

2 \_\dx J \dx ^ 

with the quantization condition (5.20), and apply the bosonization/fermionization for- 
mulas on H. In the Schrodinger picture, all operators are of course regularized by point 
splitting at equal time. Fermions can be defined by using any scalar fields X and X 
satisfying the same quantization condition as and (5.20) by 



27r V 27r 





For our purposes, we need 



The current j'^ — ip^^ip is 



X^^<P, X = (5.34) 



/ = i(V'+'0+ - = — ^ ' = -«(V'+'0+ + '0-'0-) = — (5.35) 

TT ax Qes ax 

and this is not a Lorentz vector. The four-fermions interaction is then defined by 

iPri^iP^^^iP = qoiff - qiij'f , (5.36) 

where qo and qi arc parameters that wc will adjust later to make the interaction a 
Lorentz scalar. It is straightforward to check that the fermion hamiltonian 



ax ax ji \ / n \ / 

G 



+^ {qoiff - QiW?) (5.37) 
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is equivalent to the bosonic hamiltonian 



TT 



TT 



d(f)^ 
dx 



+Fe cos{2ges(j)) + Fm cos{27r(j)/ Qes) ■ (5.38) 



The four-fermions interaction will be a Lorentz scalar if the G 7^ theory is Lorentz 
invariant with the same speed of light c = 1 as the free G — fermionic theory. This 
in turn is equivalent to the condition that i^bos coincide with H (5.32), which happens 
for 

qi = —qo (5.39) 

71 

and 

^ = 1^. (5.40) 

Equation (5.39) shows that we have one free parameter go in defining the scalar four- 
fermions interaction, and also that this interaction depends implicitly on ^fgff through qi. 
Equation (5.40) generalizes the standard Coleman's formula (5.26) which corresponds 
to the choice qo = 1. Of course the physics does not depend on go- We can now deduce 
the correct duality transformations of the four-fermions term. Equations (5.36), (5.39) 
and (5.40) implies 

G^T'^^^T.^ = [iff - ^ iff). (5.41) 

Equations (5.35), (5.27) and (5.30) show that 

3'^3\ (5.42) 

and thus by using (5.18) we finally deduce that the four-fermions term is actually 
invariant under the 5" duality. 

It would be interesting to study possible strongly coupled fixed points in theo- 
ries of the same kind as (5.17), but including other electric and magnetic charges, 
multi-components fields, and 9 angles [58]. The fact that the Argyres-Douglas CFTs 
appearing in four dimensional supersymmetric gauge theories admit an ADE classi- 
fication [59] suggests that all the minimal unitary CFTs, that also follow the ADE 
pattern [60], could be obtained in this way. A very simple example is the case where 
the electric charge = 2 of our model is replaced by an electric charge = 3. This 
yields a strongly coupled fixed point of the Potts type. 



49 



6. The double scaling limits 



We have discussed three different descriptions of the non-trivial low energy physics 
appearing in our model near the singularities of parameter space. The Ginzburg- 
Landau description (Sections 3.5 and 4.2) and the fermionic description (Section 5.3) 
are very simple, and can exist because of peculiarities of two-dimensional physics. The 
third description (Section 5.2), in terms of a theory of electric and magnetic charges, 
is a priori more generic, and I have emphasized the strong similarity with Argy res- 
Douglas CFTs in four dimensions in particular. We would like now to derive a fourth 
description, whose counterpart in four dimensions would be a string theory description 
of the vicinity of an Argy res- Douglas CFT (or more generally of the vicinity of generic 
non-trivial critical points that can be found in the context of gauge theories with Higgs 
fields). We will be rather brief on this interesting aspect of our model here, since 
we hope to provide more details in [14]. We are simply going to show explicitly that 
a double scaling limit in the sense of the "old" matrix models [13] can be defined 
when we approach the Ising or Ashkin- Teller singularities on 7iq, and that the double 
scaled theory coincide with the low energy theory in the vicinity of TYq. This provides a 
duality between the interacting theory describing the low energy degrees of freedom and 
a theory of randomly branched polymers, as explained in [61] and as will be reviewed 
in [14]. The discussion is similar to the case of the ordinary vector models, as reviewed 
for example by J. Zinn- Justin in [38]. 

There are three different types of critical behaviour occuring in our model (Ising, 
Ashkin- Teller and compactified boson), and thus a priori three different types of double 
scaling limits can be considered. Let us first discuss the vicinity of the Ising critical 
point. As discussed in Section 3.5, when we approach the critical hypersurface TCq, 
Sv — > 0, the standard 1/A^ expansion in plagued by IR divergences. This means that in 
the large N expansion of physical quantities, like (3.1), the coefficients Wi will diverge 
as 5v — > 0. Proving that a double scaling limit can be defined when Sv ^ amounts 
to proving that these divergences have a specific form for any / such that they can be 
compensated for by taking N ^ oo and Sv ^ in a correlated way. The double scaled 
theory then gets contributions from all orders in 1/A^. This means in particular that 
we must take into account the renormalization of the original non-linear a model to 
all orders. Remarkably enough, this can be done very simply. The crucial point is of 
course that the divergences in Wi are due to IR effects. 
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We have shown in Section 3.5 that at momenta p <^ A, the physics is correctly 
described by the action (3.48). By using the rescaled variables (3.51), we deduce that 
physical quantities can be correctly described by the action 

Six] = J d'x' (i ff^xd'^x -^X' + ^Vx') (6.1) 

for momenta 

p' < AVn. (6.2) 
If we now take the limit N oo, the effective UV cut-off of (6.1), which is of order 

Ao,eff ~ AVn, (6.3) 

goes to infinity. The action (6.1) thus needs to be renormalized accordingly, in order 
to get a finite limit. This is done with the help of the standard normal ordering, 

X =:x ^+6^1^^X +3(^ln^j • (6.4) 

By replacing into (6.1), we see that, up to some trivial factors, we obtain a finite double 
scaled partition function in the limit 

Sv 

N ^oo, Sv^Q, N y - 3]nN ^ constant, (6.5) 



as was to be shown. In momentum dependent quantities, p' = pyN must be held fixed. 
The same reasoning show that the correct double scaling limit in the vicinity of the 
Ashkin-Teller critical point (4.10) is defined by 

N ^ oo , Sv ^ , vn-1 ^ , 

Sv Sv V 

- AlnN ^ constant, N ^ - 4 In = constant'. (6.6) 

V V 



The last case that wc have to consider is the vicinity of the compactified boson CFT, 
for example vn-i — * with Hila^ > A^^^"^^. We cannot define a double scaling limit 
in this case. Indeed, in the N ^ oo limit, the physics is well described by the action 
(4.2), since g^Q ~ 1/A^ is then very small. We see that the 1/A^ corrections are simply 
corrections to the radius of the compactified boson, which correspond to a marginal 
operator. This is very different to the case of the Ising or Ashkin-Teller CFTs, and does 
not produce the wild IR divergences that are necessary to be able to define a double 
scaling limit. 
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A. General solitary waves solutions 

As discussed in Section 2.3.3, the most general static, finite energy solutions to the field 
equations (2.32) satisfy 

' {a-h,)^, (A.l) 



with 



5,: 



hm $2 = (A.2) 



when Hn — maxi<i<jv /ii, which we shall assume in the following, a is a Lagrange 
multiplier implementing the constraint 

^ 1 



1=1 



The set of equations (A.l, A. 3) can be interpreted as the Newton equations for the 
motion of a particle constrained to move on the sphere 5^"^ of radius 1/g m. the 
quadratic potential 

1 ^ 

um^+-Y,h,^i^-vm. (A.4) 

The role of the time in the mechanical problem is played by x, and V($) = — f/($) 
is the potential of the corresponding field theory. That such a mechanical analogy is 
possible for solitons in two dimensions is of course banal (see e.g. [42]). The peculiarity 
of the soliton problem is that we are interested only in the motions satisfying (A.2). 

We show in this Appendix that all the solitonic solutions to our model can be found, 
and are expressed in terms of elementary functions only. This generalizes the similar 
result for the sine-Gordon theory. This is possible because the classical mechanical 
problem (A.l, A. 3) is integrable in the Liouville sense. It was first studied by C. 
Neumann in the case = 3, who shows that one can separate the variables in the 
Hamilton- Jacobi equation [36], and discussed much later in the general case [62]. In 
the following, we present a self-contained elementary analysis in the case where all the 
hiS are distinct. The more symmetric cases where some of the /ijS would coincide can 
be easily obtained by taking suitable limits, but we will not spell out these details. In 
A.l we explain how to solve the equations of motion in general, and then we specialize 
to the soliton problem. 
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A.l General analysis 

The starting point are the formulas for the constants of the motion [62] , 



where = d^i/dx. Only — 1 of the XjS are independent because of the identity 

N 



^ 1 



i=l 



It is straightforward to check that dXi/dx = by using (A.l) and (A. 3). These A^ — 1 
independent integrals of the motion actually form an involutive system with respect 
to the Poisson bracket (by identifying $^ with the conjugate momentum), which prove 
Liouville integrability The energy is simply expressed in terms of the XjS, 



TV n JV 



i=i i=i 

We are thus left with N — 1 independent first order differential equations, instead of 
the second order equations (A.l). In order to solve these equations, it is very useful to 
rewrite the definition of the 2jS in a more compact, elegant way. For this purpose, we 
introduce an arbitrary complex parameter z and define 

AT 

q{x, y;z) = ^ , q{x; z) ^ q(x, x; z), (A.8) 

i=i ^ 



N 



The formulas (A. 5) are equivalent to the equation 



J($, z) = g($; z) (g($'; z) + \) - g($, zf, (A.IO) 



9 

and the differential equations we want to solve can be written 

N 



X(*.*';.).g-^.JjM. (A.n, 
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where 

TV 

a{z) = l[{z-hi), (A.12) 

i=l 

and b{z) is defined by (A. 11). The arbitrary parameter z is then chosen to make the 
equations (A. 11) as simple as possible. (A. 10) suggests that a good choice is 

q{^;z = n) = 0. (A.13) 

This equation has generically N — 1 solutions z = /ii, . . . , /ijv-i that can be expressed 
in terms of the N — 1 independent coordinates on the sphere, and vice-versa. The 
/liS are usually called elliptic coordinates. Explicitly, by introducing 

N-l 

m{z;ii) = Yl{z- iii), (A.14) 



we have 



and 



q{^-z) = ^^j^ A.15 



N-l 



n (^^-^^o 

We now make the change of variable $ i— > in (A.ll). By taking the derivative of 
(A.13) with respect to x we get 

1 A $2 

and by taking the derivative of (A.15) with respect to z we get 

^(/i,-/l,)2 ^2 a{z = iii) 
Prom these relations and (A. 10, A. 11, A.13) we finally obtain 



= 2.. V (A.19) 
n ~ 



54 



with = ±1. We can put (A. 19) in a more suggestive form by using the 
Lemma: if P{z) = YYi=i{^ ~ is an arbitrary polynomial, then 



This is proven by calculating 

I az_ z 

Jc Pi^) ^ ' 

with the help of the residue theorem and by taking C to be a circle centered at ^ = 
of radius R oo. 

The lemma implies that (A. 19) is equivalent to 

A^-l N-l-j , 

E o r f W ^ = ^^'^^ 1<J<N-1, 6, = ±1. (A.22) 

The general solution of the equations of motion are then given by elliptic integrals 
associated with the genus N — 1 hyperelliptic curve 

= -Aa{x)b{x). (A.23) 

However, as already stressed, we are not interested with the most general solutions, but 
with the ones corresponding to the solitons of our field theory model. We now focus 
on this special case. 

A. 2 The solitons for general N 
A. 2.1 Bcisic equations 

The boundary conditions (A. 2) were imposed in order to have a finite mass M for the 
solitons. M is simply given by the action of the mechanical problem, 

dx(Y,<^f + v,^A, (A.24) 



where Vi = — hi as usual. For M to be finite, it is necessary to have 



li = ^> (A.25) 



that is to say. 



7V-1 



b{z) = l[{z-h,). (A.26) 



i=l 
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We will see below that (A. 25) is actually sufficient to have a finite M. The conditions 
(A. 25) thus characterize the particular motions that can be interpreted as solitary waves 
in our original theory. 

With this particular choice for b{z), (A. 22) yields 

iV-l N-l-j 7 

E^^ nT -53,idx, l<j<N-l, Q = ±l. (A.27) 

'"^ I^Jh^^^^i JJ^ {iii - hj) 

We introduce the new coordinates 



Hi = Q V^'N - IJ'i, (A. 28) 

and we finally get our basic equation that can be written in two equivalent forms 
corresponding respectively to (A. 19) and (A. 22), 

iV-l 



Af-l 



n (y^-y^) 

N-1 

i=i 



dyi = -dx, l<i<N-l, (A.29) 



E^^^Vr^^ ^ = -i,,<;x, i<j<jv-i. (A,30) 

(A. 30) can be integrated immediately in terms of elementary functions and N — 1 
constants of integration, which gives the full solution of the problem. The fact that 
only elemetary functions are involved is due to the fact that the solitonic solutions 
correspond to a highly degenerate hyperelliptic curve (A. 23). Explicit formulas will be 
given for the case = 3 in the next subsection. Note that though (A. 30) is best suited 
to an explicit integration, a qualitative analysis of the solution, to which we now turn, 
is most easily performed by using (A.29). 

A. 2. 2 Description of the solutions 

We choose without loss of generality > h^-i > ■ ■ ■ > hi. Equation (A. 16) shows 
that $^ > implies that each interval [hi, hi+i] contains one and only one of the fijS, 
or equivalently that each interval [vj,!),-!] contains one and only one of the yjS. We 
consider a solution such that 

lim = (A.31) 
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and we conventionally choose the ordering of the coordinates such that this boundary 
condition corresponds to 

lim fji = (A. 32) 

X— >— oo 

By looking at (A. 29), we see that dyi < when x — > — oo. More generally, as long as 
Ui e [y^Vi^i, ^/Vi], which is true at least up to some finite value of x because one can 
have only one y in each of these intervals, a simple recursive reasoning using (A. 29) 
shows that dyi < for all i. Two yf can eventually swap their intervals. Equation 
(A. 29) shows that it is possible to have yf — vj at some finite x as long as some other 
yl = Vj at the same x, the pole in the denominator being then compensated by a zero 
in the numerator. With our particular boundary conditions (A. 32), the only consistent 
way this can happen is by having first yN-i ~^Jvn~i and yN~2 — ^ ^Jvn-i (step 1), 
then yN-i — ^ and yN-s — ^ ^/VN-2 (step 2), and so on up to step N — 2 when 

— > and yx \fv2- This is followed by — 3 additional steps during which 

y\ -Jvl and yN-2 -^/vl (step A^ - 1), |/i ^ ^Jvl and y^-z -y/vl (step A^), and 
so on up to step 2 A" — 5 when yi — > ^Jvn-i and y2 —^/vN-l■ At each intermediate 
step it can be checked using (A. 30) that dy^ < for all i, so the y^s are monotonic 
decreasing functions oi x. At the end of the day we have 

lim yi = -^VN-i- (A. 33) 

What we have been describing is the generic solution. By restricting oneself to motions 
for which some of the $jS are equal to zero, one can obtain other solutions, correspond- 
ing to the generic solution at a lower value of A^. 

To get a better physical understanding of the nature of the solution, let us consider 
the coordinates $i. First, let us note that the transition between the different steps 
correspond to the successive crossing of the hyperplanes $7v-i = 0, ^n-2 = 0, . . . , $2 = 
0, $3 = 0, . . . , ^N-i — 0. This shows that the trajectory lies entirely on the half sphere 
$1 > or $1 < 0, the two cases being related by the symmetry ^2(1)- Second, and 
more importantly, as x increases we have successively yN-i — 0, yN-2 — 0, and so on 
up to yi = 0. This shows that 

^% = \ (A.34) 

has N — 1 zeros when x goes from —00 to +00. Prom this we can conclude that when 
A^ is odd the solution is topologically trivial, while when N is even, it belongs to the 
topologically non-trivial sector of the theory. Equivalently, we have 

x—*+oo g 
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The fact that we cross the equator — N — 1 times suggests that the solution actu- 
ally corresponds to the succession oi N — 1 kinks and anti-kinks. The N — 1 parameters 
(constants of integration) would be in this interpretation the global center of mass and 
the relative separations of the individual kinks. To prove that this interpretation is 
indeed correct, we are going to show that the total mass (A. 24) of our solution is the 
sum of the masses of the individual sine-Gordon solitons of the model (2.56), 

2 ^"^ 
^ 1=1 

Note that M does not depend of the N —1 parameters of the solution, which is a direct 
consequence of the fact that M coincide with the action of the mechanical problem, and 
that the boundary conditions in the action are independent of those parameters. We 
have thus obtained a static solution consisting of a succession of sine-Gordon kinks and 
anti-kinks! As pointed out before, we actually have solutions for any number k < N — 1 
kinks and anti-kinks. Note that all the kinks or anti-kinks appearing in the solution 
have different masses. A way to prove (A. 36) is to choose the constants of integration 
in some limit that simplify the general solution (kinks and anti-kinks in a limit of very 
large separation). We will present instead a direct derivation which turns out to be 
simpler. 

Our starting point is 



„+oo N-l 

M= / dxJ2^i^l (A-37) 

-^-"o i=l 

which is derived from (A. 24, A. 25, A. 7). We then use the following algebraic identities: 

7V-1 Af-1 

i=l ^ i=l 

N-l 



Y.i^^-y^)-Y.^i (A-39) 

n ^vl-y^) 



N-l N-l 



1=1 



j=l,j^i 

Equation (A.38) is proven by writing (A. 15) as 

^ m{z; ii) = a{z)q{<^] z) (A.40) 

and identifying on both sides the coefficient of z^~'^. Equation (A. 39) is proven by 
using the fact that the right hand side can be viewed as a rational function of yl of 
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degree d < 1, and that it is actually a linear function of yf since it has a finite limit 
when yl ^ yf , ior all i — 2, . . . N — 1. Taking into account the permutation symmetry 
amongst the y^s, this leaves us with two unknown coefficients which are determined 
easily by looking at the y\ ^ oo limit up to terms 0{l/yf). By using successively 
(A.37), (A.38), (A.39) and (A.29), we get 

2 ^'^ r 

^ = --2Y. dm. (A.41) 

^ 1=1 

By using the fact that the y^s are monotonic functions of x, and (A. 32) and (A. 33), 
wc finally obtain the desired result (A. 36). From our understanding of the solution, we 
deduce in particular that the mass density 



1=1 

will have N — 1 maxima, located at the positions of the kinks and anti-kinks, whose 
values approach the corresponding values for sine-Gordon solitons, Vi/g"^, at large sep- 
arations. 

A. 3 Formulas for the case N=3 

We give below the explicit formulas used in Section 2.3.3 to make the Figure 2. By 
integrating (A. 30) for TV = 3 we get 

(V^-yi)(V^-l/2) ^ g2V^(x-a.i) - V^)iV^-y2) ^ ^2^(x-x,) 

i^M + yl){^/vl + y2) ' {^/v^ + yl){^/v^ + y2) 

When I ^2 — 3:^1 1 >> 1/ ■s/v^, we have two well separated kinks centered aX x — xi and 
X — X2- Solving explicitly for yi and y2 amounts to solving a degree 2 polynomial 
equation. The relevant root is picked up by using (A. 32); the resulting formulas are 
very complicated and will not be listed here. The spherical coordinates (6*, (p) are then 
given by 



arccos ^^^^ i = signfyi — Jv2j arctan W — ^ — ^^^^^^ — (A. 44) 
y/viv^ V V2 {vi - y(){vi - y|) 
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B. The functional s[f] 

In this Appendix, we study the functional s[f] defined by 

s[f]='-tTln{-d,d^-f + te) + j^ln^ Jd'xf{x), (B.l) 

or equivalently in the euchdean by 

= ^ trln(-a2 + /)--! In ^ [d'xf{x). (B.2) 
2 An ji J 

This functional is ubiquitous in the study of the large limit of our model. s[f\ 
depends on an arbitrary renormalization scale 

All UV divergent integrals are regulated by a momentum cutoff \p\ < Aq. In finite 
quantities, we will take the limit Aq — > oo. 

B.l s\f] for constant / 

For constant /, s[f] can be written in terms of a potential v[f) as 

s[/ = cst] = Jd^xv{f). (B.3) 
When / is a positive constant, the euclidean formula (B.2) yields 

which is, up to /-independent terms. 

The correct formula for / < is obtained by using Feynman's ie prescription (B.l), 

^(/) = -^(ln^-^7r), for/<0. (B.6) 
Equation (B.6) is the analytic continuation of (B.5) by going to Im/ < 0. 



B.2 s[f] for arbitrary / 

When / is arbitrary, one can expand s[f] in terms of ordinary one-loop Feynman 
diagrams. Introducing an arbitrary mass scale m and defining (f){x) = f{x) — m^, we 
have 

oo 

s[f] = Y,Sn[f;m% (B.7) 



n=0 
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where 



1 /■ " 

J i=l 

1=1 ^ ' i=l 



with 



(2x)V(X:p.) 5<">(p„...,p„;m2) = 

/n 
JJ(i2a;ie^^?=i*'^-^^s(")(a;i, . . . ,a;„;m2). (B.IO) 



The euchdean formula for s*^"^ is 



i=0 j=l 

Exphcitly, the hnear term is given by 

1 , 



(B.ll) 



5« = -!-ln^ (B.12) 



or equivalently 

si[f;m']^^\n^ f£x<t>{x). (B.13) 
In the euchdean regime > 0, the quadratic term is given by 

s^^) (p, -p; m^) = s^^) (p^ m^) (B.14) 

1 1 ^ v/l + 4mVp^ + l 

= — , = m ^ , , for p > 0. 

47r p2yi + 4mVp2 ^1 + 4m7p2 - 1 

At low momentum we can use the expansion 
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or equivalently 

''^^''"'^^"leW /^'"''^'^"'^"geW Jd'xcf>{x)cf4>{x) + o{&'). (B.ie) 

For physical momenta below the pair production threshold, < —p^ < 4m^, the correct 
analytic continuation is given by 

s^'^\p'^] rn^) — -^^^^=^^^= arctan — -^^^^=^^^= , for < —p^ < 4m^. 

27r p2^-i- ^rin?/p'^ y-l - 'irin?/p^ 

(B.17) 

Finally, for —p^ > 4m^, which is the region relevant for two-particle scattering, we have 

~(2)t 2 2^ 1 1 A 1 + yi + 4mVp' , 2,2 

' [p ,m ) = = m ^-^^^=^^= — ZTT , for — p > Am . 

47r pVl + 4mVp2 \^ 1 - yi + 4mVp2 J 

(B.18) 

In the latter case, it is convenient to introduce the rapidity parameter ^ > such that 

-p^ = 4m^ cosh^(^/2) (B.19) 

and in terms of which 

Svrm^ smht^ 

Let us quote to finish the asymptotic behaviour of s*-^-* at large euclidean p^, which is 
useful to investigate the UV properties of the 1/N expansion, 

= -7^ (in 4 + ^ (l - 1- 4) + >n 4)) ■ (B-^D 

47rp V m p \ m / \ p m / / 

C. The 1/N corrections to the critical hypersurface 

We discuss in this Appendix the corrections to the equation 

N-l 
i=l 

for the hypersurface of singularities Tiq derived in Section 3.4. The main goal of this 
calculation is to understand the role of the infrared divergences, which generically 
plague the 1/N expansion near Tiq. As the very existence of Tiq was derived within the 
large expansion, the consistency of the analysis requires that the equation for Tiq 
itself gets only small corrections when we go beyond the leading approximation. We 
show below that the first corrections are of order A^~^ In A^, larger than the naive A^~^, 
but nevertheless smaller that the leading A"° term. The calculation also illustrates nicely 
how renormalization works within the 1/N expansion, consistently with the discussion 
of Section 3.1. 
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C.l The effective potential 



The starting point is the effective action 



S'efr[ci;, $Ar] = (A^ - 1) Ses[a, (f] 



(C.2) 



where 



SeS 



[a, (p\ — J d^x 



— (C.3) 

1=1 



We were careful in keeping subleading terms in (C.3), since we are now willing to take 
them into account. We have defined 



(p = ^n/VN - 1 



(C.4) 



and A by the equation 



1 Ao 
g'o 27r A ■ 



(C.5) 



To order we have A = A, but to order N~^, A must pick an infinite, Ao-dependent 
contribution, coming from the 1/N corrections to the P function. Equations (3.4), (3.7) 
and (3.8) of Section 3.1 gives the form of this contribution, as well as the renormalization 
of hifi — hiZi^s,2)/ Z: 



= A^ 



l + lln^-llnln^ + ^ 
N N AN 



(C.6) 
(C.7) 



where Ci and C2 are finite constants to be determined by some renormalization condi- 
tions. 

The 1/N corrections correspond to the one- loop diagrams derived from the non- 
local action Seff- In particular, the effective potential Ves to order N~^ is such that 



A 

N-l 



i-^ X:(«-M 1^^^) +^trlnj[a,^]. (C.8) 
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J is the functional hessian of Seff , for which 



(C.9) 



where we have defined 

7V-1 

~s{j?- a,h,,..., hN-i) = 5^ ^(^Hp'; « - ^0- (c.io) 

i=l 

Prom (C.8, C.9) we get the basic formulas, valid up to terms of order l/N'^, 

- {a-hM,o)v-^--:, dp-—^ —~— ^ r'(C.ll) 



'"'^ NAtt'J {p' + a-hN)S{p';a,hi,...,hN-i)-ip' 

dves ^ ^ 1 , a - hj^ 1 C ^ S + {p'^ + a - hN)daS 

da ~ 2 87r(7V-l)^" ^ Stt^TV i ^ (p2 + ^ _ _ ^2 " ^ 

C.2 The equation for Hq 

The naive strategy to get the equation for Hq up to order 1/N would be to solve 

dv^s dv^s 



dip da 



(C.13) 



by substituting a and (/? in the terms of order by their respective values at order 
iV°, that is a = hj^ and (p = (see Section 3). However, this yields IR divergent 
integrals, and thus is not correct. Instead, we will identify the terms responsible for the 
IR divergences, and set a — and (p — in the other terms only. The IR divergences 
can be analysed easily by using 

S{p^; a^hN,hi,..., hN-i) = + O(p'), (C.14) 

OTTV 

1 1 

dJip'; a^hN,h,,..., hN-i) = g^.^ _ E ^ + ^(^'')' 

where V is defined in (3.49). The term containing daS in dv^s/da is not IR divergent, 
while the other term can be written 
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8^ 



S{p'^;a,hi,. . . ,hN_i) 



+ 



(p2 + a - hN)S{p^; a,hi,..., Hn-i) - 

Sjp"^ = 0;a,hi, . . . , hN-i) 

(p2 + a - hN)S{p^ = 0;a,hi,. . ., Hn-i) - (p^ 

f (Pp_ = 0;a,hi, . . . , hN-i) 

J 87r2 (p2 + ^ _ hN)S{p^ = 0; a, /ii, . . . , /ijv-i) - V'^ 



(C.16) 



The first term is now IR convergent, and goes actually to zero when o; — > /ijv and 
</? — > 0, while the second term is easily calculated. The condition dves/da = is thus 
equivalent to 

^ ^ _ 1 a-/l^,o 1 Ag Ia{vi,...,VN-l) 

- ' V - 1) ' " -r „ ni , „ _ „ 



2 87r(A^-l) ^ A2 SttA^ q; - /ijv + SttT/c/p^ 



SttA^ 



where 



Iaivi,...,VN-l) = 3 / 



daS{p'^] a = hN,hi, . . . , hN-i) 
S{p^;a = hN,hi,. . . ,hN-i) 
We treat dv^s/dip is a similar way. We write the 1/N contribution as 



1 



(Pp 

4^ 



1 



(p2 + q; - hN)S{p^] a,hi,..., hN-i) - (fi'^ 

1 

(p2 + a - hN)S{p^ = 0;a,hi,..., /ijv-i) - 
1 



47r2 (p2 + a - /iiv)5(p2 = 0; a, /ii, . . . , /ijv-i) - 



(C.17) 
(C.18) 



(C.19) 



check that the first term is now IR convergent while the second term is easy to compute, 
and conclude that the condition dves /d(fi — is equivalent to 



V / 

= a- Vo + T7 21n , °^ , 



- I^{vi,...,vn~i)^ , 



(C.20) 



where 



, ( _J_ ffp ( 1 



= 0; a = hN, h, . . . , /ijv-i) 



(C.21) 
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We have chosen to approach Hq from weak couphng, and thus we have discarded the 
solution (fi — to the equation dves/d(fi — 0. Tiq is precisely the locus in parameter 
space where (p — becomes the true minimum of the effective potential for (p, see 
Section 3.4. 

The equations (C.17) and (C.20) arc our new starting point. They are badly UV 
divergent, as discussed further in the next subsection, and the renormalizations (C.6) 
and (C.7) arc not enough to make them finite. This is not surprising: as explained in 
Section 2.1, formulas containing the field a need an infinite number of counterterms. 
To obtain a finite, physically sensible, formula, we must eliminate a using (C.17), and 
then check that the result is finite. To make the presentation as clear as possible, let 
us introduce 

a^a- hNfi- (C.22) 

We know from Section 3.4 that a — on Hq to order (this is a direct consequence 
of (C.20). We will thus always set a = in 1/N corrections, as long as this does 
not lead to an IR divergence. Moreover, to leading order, (C.17) shows that a — 
&No{vi,o, . . . , vn-i,o; A) with 

N-l 

y In "^\+ ^ 0. (C.23) 
1=1 

In order to solve (C.17), we then substitute a — oino-\-0{1/N) and use (C.23) to obtain 

« = «^.a + ^ (/a - In ^) . (C.24) 

Finally, using (C.20), we obtain the ct-independent equation 

Q!ivo (vi,o, . . . , VN-ifi] A) = 1^/^-4 + 3 In j . (C.25) 

This is the equation for TYq at order 1/N. A highly non-trivial check is that it is 
actually Ao-independent. A crucial ingredient for this to be possible is that only the 

combination — /„ appears, so that most of the UV divergences in and la cancel 
each other (sec (C.41) for an explicit formula). Moreover, by using (C.6) and (C.7), it 
is straightforward to show that up to terms of order l/N^, 

OiNo{vifi, . . .,vn-i,o; A) = aNo{vi, . . . , wat-i; A) + — (in-^ - 41nln-^ + ci - C2 j . 

(C.26) 

This turns out to be exactly what is required to cancel the divergences in (C.25). In 
terms of finite quantities only, the equation for Tiq can then be written 

V f a NO 



a 1^0 



= -(/-41n2 + 31n^ + C2-Ci). (C.27) 
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The finite integral I{vi, . . . vn-i) is defined by (C.42). It remains to enforce a renormal- 
ization condition in order to determine C2 — ci. The study of this seemingly secondary 
point is actually important, because the correction in (C.27) are small only if we can 
choose C2 — Ci to be much smaller than A^. In particular, the naive guess that C2 — Ci 
can be an arbitrary finite constant independent of N is not correct. Our renormaliza- 
tion condition will be written at Vi — ■ ■ ■ — vn-i — v, values for which / = and 
a^o — —V + A^. We impose that 



A' 



(C.28) 



A being an A?"- independent, strictly positive, constant. It is impossible to choose A = 
0, due to the IR divergences, very much like it is impossible to use renormalization 
conditions at zero momentum in massless theories. Equation (C.28) implies 

C2-C1 = 31nA^ + A-31nA + 41n2, (C.29) 

and the final equation for Tiq is 

^(sinAT + J + A + Sln^). (C.30) 
We recall that a^o is the unique solution to 



V 



N-1 



^In T^ = 0, (C.31) 



1=1 



A2 



/ is defined by (C.42), and A > is a renormalization constant independent of A^. We 
thus see that the IR instabihty is responsible for a correction of order (In N)/N, larger 
than the expected 1/A^ but much smaller than the leading N° term, as required. 

C.3 Formulas for and — 

In view of (B.14), it is natural to introduce the variables 

..(p^) ^ (0.32) 

This set of variables is particularly well suited to study the integrals and when 
the ViS are equal or nearly equal. The UV behaviour is also easily studied in this 
representation. One has 
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with 



so that 



^ 1 / , / 1 + ^^'"^ ~ 

Vi {Xi + 1)2 

9ij{xk)gji{xi) = gjk{xk)- 



{CM) 



(C.35) 



We have the exphcit formulas 



f. = -E/ 



Xi-lf , 2 (x,-l 
N-i ^^.(a2\ —77-, — -\ I r I mx,; 



EXj[Xi) 1 , , 



(C.36) 



^1 /-,(Ag) 



i=l 



x^-l 



JV-1 



E 



^ + 1 



In Xj (xi) 



Xi + 1 



N - 1 Xi{xi - 1) 



(C.37) 



/(p and Jq, have separately untamable UV divergencies. These divergencies can be 



studied for example for vi 
I^{vi = ■ ■ 



vn-i. In that case one has 



with 



A2 A2 
VN_i = v) = 21n^-21nln^ + 27-ha;(A^), 

V V 



(C.38) 



(C.39) 



and li is the logarithmic integral, whose expansion is 
\ix = P 



dx , , (Inx)*^ „ 

= 7 + minx + > — TT-; — for x > 1. 



Inx 



k=l 



k\k 



(C.40) 



On the other hand, — la has simple UV divergencies that are renormalizable with a 
finite number of counterterms. The general formula is 



/^-/a = /-41nln^ + 41n^_41n2 + ^ 5]ln 



(C.41) 



i=l 
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where / is a finite integral (which vanishes when vi 



vn-i), 



I{Vi, . . .,Vn-i) = X] / d^i 



N — 1 Xilnxi 

4(x, - 1) 
Xi{Xi + 1) 



4 1 



N - 1 Xi{xi + 1) 



2 Xi-1 



1=1 




lna;i 



■ (C.42) 



On the critical hypersurface TCq, the last term in (C.41) is of order 1/N and can be 
neglected. 
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